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ABSTRACT 

The abundances of iron and oxygen are homogeneously determined in a sample of 523 nearby (d < 150 pc) FGK disk and halo stars 
with metallicities in the range -1.5 < [Fe/H] < 0.5. Iron abundances were obtained from an LTE analysis of a large set of Fei and 
Fe ii lines with reliable atomic data. Oxygen abundances were inferred from a restricted non-LTE analysis of the 777 nm O i triplet. 
We adopted the infrared flux method temperature scale and surface gravities based on Hipparcos trigonometric parallaxes. Within 
this framework, the ionization balance of iron lines is not satisfied: the mean abundances from the Fe i lines are systematically lower 
by 0.06 dex than those from the Fen lines for dwarf stars of T e s > 5500 K and [Fe/H] < 0.0, and giant stars of all temperatures 
and metallicities covered by our sample. The discrepancy worsens for cooler and metal-rich main-sequence stars. We use the stellar 
kinematics to compute the probabilities of our sample stars to be members of the thin disk, thick disk, or halo of the Galaxy. We find 
that the majority of the kinematically-selected thick-disk stars show larger [O/Fe] ratios compared to thin-disk stars while the rest show 
thin-disk abundances, which suggests that the latter are thin-disk members with unusual (hotter) kinematics. A close examination of 
this pattern for disk stars with ambiguous probabilities shows that an intermediate population with properties between those of the thin 
and thick disks does not exist, at least in the solar neighborhood. Excluding the stars with unusual kinematics, we find that thick-disk 
stars show slowly decreasing [O/Fe] ratios from about 0.5 to 0.4 in the -0.8 < [Fe/H] < -0.3 range. Using a simple model for the 
chemical evolution of the thick disk we show that this trend results directly from the metallicity dependence of the Type II supernova 
yields. At [Fe/H] > -0.3, we find no obvious indication of a sudden decrease (i.e., a 'knee') in the [O/Fe] vs. [Fe/H] pattern of 
thick-disk stars that would connect the thick and thin disk trends at a high metallicity. We conclude that Type la supernovae (SN la) 
did not contribute significantly to the chemical enrichment of the thick disk. In the -0.8 < [Fe/H] < +0.3 range, thin-disk stars show 
decreasing [O/Fe] ratios from about 0.4 to 0.0 that require a SN la contribution. The implications of these results for studies of the 
formation and evolution of the Galactic disk are discussed. 
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1. Introduction 

Nearby stars are ideal targets for chemical abundance studies due 
to the availability of accurate trigonometric parallaxes, a negligi- 
ble reddening, and their brightness, which allows us to estimate 
the stellar atmospheric parameters as well as acquire high reso- 
lution, high signal-to-noise spectra of these objects. 

With the Sun practically in the Galactic plane, the solar 
neighborhood contains mainly thin-disk objects, with a still sig- 
nificant fraction of thick-disk stars and only a few halo mem- 
bers. Distinctive patterns in the chemical composition of cer- 
tain groups of stars contain valuable information for studies of 
Galaxy formation, evolution, the chemical enrichment history of 
the interstellar medium, and the nature of the nucleosynthesis 
sites. This paper deals with the determination of iron and oxy- 
gen abundances in nearby stars. 

* Tables 4 to 6 are only available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+A/ 
** Partially based on observations obtained with the Hobby- 
Eberly Telescope, which is a joint project of the University 
of Texas at Austin, the Pennsylvania State University, Stanford 
University, Ludwig-Maximilians-Universitat Miinchen, and Georg- 
August-Universitat Gottingen; and data from the UVES Paranal 
Observatory Project (ESO DDT Program ID 266.D-5655). 



Due to the large number of lines and the availability of reli- 
able atomic data, iron is an excellent reference element to assess 
the overall chemical composition of FGK stars. In fact, the term 
"metallicity," denoted as [Fe/H], usually refers to the iron abun- 
dance on lyQOxy gen is also one of the most important elements 
in chemical abundance studies. It is the third more abundant el- 
ement in the Universe, after hydrogen and helium, and its abun- 
dance has a relatively important effect on the internal structure 
and evolution of stars. 

Unfortunately, the determination of stellar abundances is not 
straightforward. Since they are the basis upon which many as- 
trophysical studies rest, it is important to examine carefully the 
methods of abundance determination and possible errors associ- 
ated with them. Oxygen is a particularly difficult element, with 
only a few lines available to analyze and with all of them sub- 
ject to different model complications and/or observational dif- 
ficulties. For example, in main-sequence FGK stars, the forbid- 
den [Oi] lines at 630 nm and 636 nm are both very weak and 
the 630 nm line is blended with a Ni line; the UV and IR OH 
lines are both very difficult to observe, the former due to atmo- 
spheric extinction and low CCD sensitivity and the latter be- 



1 Throughout this paper, we use the standard definitions [A'/y] = 
\og(N x /NY)-log(Nx/N Y )e and A x = log(N x /N H )+ 12, where N x is the 
number density of the element X. 
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cause they are weak lines. The UV OH lines are also strongly 
affected by surface inhomogeneities (Asplund & Garcia Perez 
2001). Finally, the Oi lines at 777 nm, which we use as our 
oxygen abundance indicators, even though easy to observe, are 
severely affected by departures from LTE (e.g., Kiselman 1993). 

The global structure of the Galaxy has been revealed by anal- 
yses of star counts. In most of these studies, the star count data 
are well approximated by two double-exponential disks, for the 
thin and thick disks, and two spheroids, for the bulge and halo 
(e.g., Chen et al. 2001, Siegel et al. 2002, Larsen & Humphreys 
2003, Robin et al. 2003, Juric et al. 2005)0 Using star counts, 
Gilmore & Reid (1983) were the first to provide compelling ev- 
idence for the existence of the Galactic thick disk. They sug- 
gested vertical scale lengths of 300 pc for the thin disk and 
1450 pc for the thick disk, as well as a 2% relative density of 
thick-disk stars in the solar neighborhood. More recent studies 
have confirmed the existence of the two disks in star counts, al- 
though the values of the scale lengths (both vertical and radial) 
as well as the local relative density of thick-disk stars are still un- 
der debate (Chen et al. 2001, Siegel et al. 2002, Cabrera-Lavers 
et al. 2005). Values reported in the literature for the vertical scale 
length range from 100 to 400 pc for the thin disk and from 600 
to 1500 pc for the thick disk. The radial scale length of the thick 
disk (~ 3 kpc) also seems to be larger than that of the thin disk 
(~ 2 kpc), and estimates of the local relative density of thick 
disk stars range from 2 to 15%. Thick disks have also been ob- 
served in many edge-on spiral galaxies, where bright thin disks 
are seen surrounded by faint red thick disks in deep photometry 
(e.g., Dalcanton & Bernstein 2002). 

The thin- and thick-disk stars of our Galaxy have distinct 
Galactic space velocities: thick-disk stars appear to be rotating 
slower around the center of the Galaxy (mean drift velocity rel- 
ative to the Sun of about -50 km s~') and have wider velocity 
dispersions compared to thin-disk stars (mean drift velocity of 
about -10 km s _1 ; e.g., Soubiran 1993, Soubiran et al. 2003). 

Chemical abundance surveys demonstrate the metal-poor na- 
ture of the thick disk compared to the local thin disk (e.g., 
Gilmore et al. 1995, Allende Prieto et al. 2006). The [Fe/H] dis- 
tribution of the thick disk spans the range -1.0 < [Fe/H] < -0.3 
while that of the thin disk is -0.8 < [Fe/H] < +0.40 It seems 
that no thin-disk stars with [Fe/H] < -0.8 exist whereas the ex- 
istence of thick-disk stars with [Fe/H] > -0.3 is under debate. 

Fuhrmann (1998) has suggested an age gap, at about 9- 
10 Gyr, between thin- and thick-disk stars, such that all thick- 
disk stars are older than thin-disk stars. Later studies (e.g., 
Gratton et al. 2000, Bensby et al. 2005, Reddy et al. 2006, 
Allende Prieto et al. 2006) have confirmed the old nature of 
thick disk stars but the age gap is difficult to detect due to un- 
certainties and systematic errors in the stellar age determination. 
Note, for example, that in their study of a very large sample of 
stars from the third data release of the Sloan Digital Sky Survey, 
Allende Prieto et al. (2006) find thick-disk stars as young as 
8 Gyr. 



2 In general, the double-exponential disks have the form: p(R, z) = 
p a e~ ( - R ~ Re),hR e~ a/hz , where R is the galactocentric distance, z the height 
above the Galactic plane, p Q and R e the local values of p and R, h R the 
radial scale length, and h z the vertical scale length. As for the spheroids, 
there is a wide variety of density laws (see Sect. 6.1 in Siegel et al. 
2002). 

3 In this paper, the so-called metal-weak thick-disk is not studied due 
to the low number of stars with [Fe/H] < -1 in the sample; the reader 
is referred to Beers et at. (2002) for a discussion of this topic. 



Regarding oxygen abundance patterns, they have been 
suggested to be distinct for thin- and thick-disk stars (e.g., 
Prochaska et al. 2000, Bensby et al. 2004): thick-disk stars ap- 
pear to have larger [O/Fe] and [O/H] compared to thin disk stars 
at any given metallicity within the [Fe/H] range in common. 
This implies that thin- and thick-disk stars were formed from 
two different mixtures of gas, presumably at different epochs 
and/or locations. Oxygen abundance patterns and their system- 
atic differences between thin- and thick-disk stars thus contain 
important information about the formation and evolution of the 
Galactic disk. Other elements, but not all of them, also show 
distinct abundance patterns for thin- and thick- disk stars (e.g., 
Bensby et al. 2005, Reddy et al. 2006). 

This paper is distributed as follows: basic data are described 
in Sect. [2] the stellar kinematics and the criterion adopted for thin 
disk, thick disk, or halo membership are discussed in Sect. [3] the 
abundance analysis is explained in Sect.|4]with the exception of 
the non-LTE corrections to the oxygen abundances, which are 
discussed separately in Sect. [5] The oxygen vs. iron abundance 
trends are traced in Sect. [6] stellar ages are discussed in Sect. [7] 
and interpretations of the results are given in Sect. [8] 



2. Sample and Data Sources 

2.1. Spectra 

Most of the spectra used in this work were available from pre- 
vious and ongoing studies of nearby stars and were thus already 
reduced using standard methods (see the references below for 
details). In addition, we obtained spectra of 22 stars to comple- 
ment the available data. A summary of our data sources is given 
in Table Q] 

Spectra of about 100 stars within 15 pc from the Sun 
were available from the Spectroscopic Survey of Stars in the 
Solar Neighborhood (S 4 N) of Allende Prieto et al. (2004a, 
hereafter AP04a)Q Their observations were made with the 
2.7 m Telescope at McDonald Observatory and the ESO 1 .5 m 
Telescope at La Silla. Spectra of about 50 stars obtained with the 
9.2 m Hobby Eberly Telescope (HET) as part of ongoing studies 
of the closest stars, metal-rich dwarfs and giants were also avail- 
able for this work. We also used 44 FGK stars spectra from the 
online UVES-VLT library by Bagnulo et al. (2003)0 who took 
advantage of the availability of the UVES spectrograph on the 
VLT during part of twilight time to observe a significant number 
of bright stars covering most of the HR diagram as well as stars 
in two open clusters. 

The major contribution of data to this work comes from the 
Reddy et al. (2003, 2006; hereafter R03, R06) studies. About 
180 of these spectra were used in the R03 chemical abundance 
work that consists mainly of thin-disk objects, while spectra of 
about 170 stars with halo and thick-disk kinematics were kindly 
provided by B. Reddy prior to publication (R06). 

Due to the lack of metal-rich stars with thick-disk kinemat- 
ics in the samples described above, we obtained spectra of 22 
stars with these characteristics. We used the 2dcoude spectro- 
graph on the Harlan J. Smith 2.7-m Telescope at McDonald 
Observatory (Tull et al. 1995) with a resolving power of about 
60,000. This is the same instrument used for the Reddy et al. 



4 http://hebe.as.utexas.edu/s4n, 
http : / /www . astro . uu . se/~s4n 

5 http: /www. eso.org/uvespop 
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Source 


Facility/Instrument 


R 


S/N 


Spectra 


AP04a 


McD/2dcoude 


60,000 


150-600 


102 


ESO/FEROS 


45,000 


TW 


HET/HRS 


120,000 


~ 300 


49 


B03 


VLT/UVES 


80,000 


300-500 


44 


R03 


McD/2dcoude 


60,000 


-400 


178 


R06 


McD/2dcoude 


60,000 


100-200 


169 


TW 


McD/2dcoude 


60,000 


~ 250 


22 



Table 1. Sources of data (spectra only) used in this work: 
AP04a (Allende Prieto et al. 2004a), B03 (Bagnulo et al. 
2003), R03 (Reddy et al. 2003), R06 (Reddy et al. 2006), TW 
(this work). Facilities: McD (2.7 m Telescope at McDonald 
Observatory), ESO (1.5 m Telescope on La Silla), HET (Hobby- 
Eberly Telescope), VLT (Very Large Telescope). 
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Fig. 1. Distribution of our sample stars in V magnitude and dis- 
tance (d). Dash-dotted lines: P\ > 0.5 (i.e. mostly thin-disk 
stars), solid lines: P2 > 0.5 (thick-disk stars), dotted lines: 
P3 > 0.5 (halo stars). 



studies. Overscan subtraction, fiat-fielding, scattered light re- 
moval, extraction, wavelength calibration, and continuum nor- 
malization were applied to the spectra using IRAF0 

Combining all data sets we ended up with 564 spectra for 
523 stars, after discarding high rotational velocity stars, double- 
lined spectroscopic binaries, and very low S/N spectra. All spec- 
tra we used have high-resolution (R ~ 45, 000- 120, 000), mod- 
erate to high signal-to-noise ratios (S /N ~ 100 - 600), and cover 
at least most of the wavelength range required for this work 
(4500 A - 7800 A). For the stars that have spectra from more 
than one source, due to different values of the resolution, instead 
of combining the spectra, the analysis was done independently 
and the final results averaged. 

To perform a differential abundance analysis, a solar spec- 
trum was adopted from the S 4 N database. This spectrum was 
obtained by co-adding skylight spectra taken when the angle 
between the Sun and the port that allows skylight to enter the 
slit room (directed towards the zenith) was between 35 and 60 
degrees, which minimizes the distortion of the solar integrated 
spectrum due to scattering in the Earth's atmosphere (see, e.g., 
Gray et al. 2000). A comparison of this spectrum with the solar 
flux atlas of Kurucz et al. (1984), smoothed to match the reso- 
lution of the former, showed no evidence of dispersed light and 
that the observational error is instead dominated by the contin- 
uum placement (Sect. 2 in AP04a). In addition, the reflection 
spectrum of the asteroid Iris, used by R03, as well as another 
skylight spectrum from our latest observations, were used. 

2.2. Photometry 

Stellar magnitudes and colors were required to obtain surface 
gravities and effective temperatures from the location of the 
stars in the HR diagram and temperature-color calibrations, re- 
spectively (Sect. 14.3b . The General Catalogue of Photometric 
Data (Mermilliod et al. 1997) provided magnitudes and col- 
ors of our sample stars in several systems (UBV, uvby, Vilnius, 
Geneva, RI(C), and DDO). Photometric data from the Hipparcos- 
Tycho mission (ESA 1997) and the Two Micron All Sky Survey 
(2MASS, Cutri et al. 2003) were also used. 

All our sample stars are in the Hipparcos catalog and have 
parallaxes measured with uncertainties typically smaller than 
10% (the average of these parallax errors is 5.8%). The distri- 
bution of our sample stars in V magnitude and distance is given 

6 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation - http : // iraf . noao . edu 



in Fig. [TJ In this Figure, the stars have been assigned to either 
the thin disk, thick disk, or halo of the Galaxy, as described in 
Sect. 13.31 Most of our sample stars are brighter than V ~ 10 mag 
and closer than d ~ 100 pc. Note, however, that in order to col- 
lect a significant fraction of halo and thick-disk stars, a larger 
volume needs to be sampled and this is why, in general, our halo 
and thick-disk objects are fainter and more distant than the bulk 
of our thin-disk members. 

In Fig. [TJwe see that most (91%) of our stars are within 
100 pc from the Sun. Therefore, for most of them, reddening 
is not an issue. For the farthest stars, however, corrections of the 
order of 0.01 - 0.02 mag may be required if they are close to the 
Galactic plane, but they were not applied here. There is virtually 
no interstellar absorption within 75 pc from the Sun (see, e.g., 
Fig. 3 in Lallement et al. 2003). This may also be the case for 
larger distances as the shape of the 'Local Bubble' (a region of 
very low interstellar gas density, devoid of dust, around the Sun) 
is far from spherical (e.g., Leroy 1999) and stars as distant as 
150 pc may not need a reddening correction, depending on their 
line of sight. 

Using a combination of interstellar reddening maps and em- 
pirical laws, as described in Ramirez & Melendez (2005a, then- 
Sect. 3.1), we estimate a reddening correction E(B - V) a 
0.01 mag for stars with d 100 pc located in the Galactic plane. 
For d =i 125 pc in the Galactic plane, E(B - V) 0.02. If con- 
sidered, these reddening estimates of 0.01 and 0.02 will lead to 
higher effective temperatures by about 50 K and 100 K, respec- 
tively, at r e ff a 6000 K. R06 applied reddening corrections to 
stars farther than 100 pc and found the maximum E(B - V) to be 
about 0.03 (roughly equivalent to their quoted Ay ~ 0.1), which 
corresponds to a maximum increase in r e jf of about 140 K. Note 
from Fig. [TJthat our sample contains a comparable number of 
thin- and thick-disk stars beyond 100 pc which implies that the 
errors introduced by neglecting reddening affect in a similar way 
our thin- and thick-disk subsamples. 

3. Kinematics 

3.1. Radial velocities 

For most (94. 1 %) of our sample stars, we adopted radial veloc- 
ities (V r ) from the literature. A careful inspection of the values 
given by different authors allowed us to discard a few published 
results and/or assign realistic error bars to the literature data. 

For 75.9% of the sample we found V, values in the 
Nordstrom et al. (2004) catalog and adopted them. Then, the 
radial velocitites given by AP04a and R06 were adopted for 
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another 1 1 .9% of the sample stars that are not included in the 
Nordstrom et al. catalog. Note that AP04a adopted mean val- 
ues of literature data while R06 used published results and, for 
the stars without a previous determination of V r , derived ra- 
dial velocities by measuring wavelength shifts for a large set of 
spectral lines. For 6.3% of our sample stars, which are not in- 
cluded in any of the studies cited above, we adopted the mean 
V r from the following sources, when available (we checked and 
avoided duplicity of original results among catalogs when com- 
puting the mean): Dufiot et al. (1995), Barbier-Brossat & Figon 
(2000), Malaroda et al. (2001), Latham et al. (2002), and Valenti 
& Fischer (2005). Finally, for only 31 stars (5.9% of the sam- 
ple) for which we did not find published V,- values, we estimated 
radial velocities as follows: first, the central wavelength of 10 
strong lines within the spectral range from 6000 A to 8000 A 
were measured and compared with those measured in the so- 
lar spectrum of Kurucz et al. (1984); the wavelength shifts were 
then transformed into velocities that include the Earth's motion 
around the Sun; and finally, these were put into the heliocentric 
frame by computing heliocentric velocity corrections using the 
rvcorrect task in IRAF. 

The mean error in our adopted V r values is 0.8 km s which 
reflects the accuracy of the literature data. For 91% of our sam- 
ple stars the error is smaller than 1 km s and it is smaller than 
5 km s _1 for 97.7% of them. Only 5 stars have errors in V r that 
exceed 10 km s _1 . These are stars with discrepant literature data, 
most likely single-lined binaries for which no orbit has been de- 
termined yet. 

We used the radial velocities to compute Galactic space ve- 
locities and determine the membership probabilities of our sam- 
ple stars to be thin-disk, thick-disk, or halo members. The accu- 
racy of our adopted V, values is sufficient for the present pur- 
poses (see below). 



3.2. Galactic space velocities 

Using Hipparcos data (coordinates, parallaxes, and proper mo- 
tions) and the radial velocities obtained as described above we 
calculated the Galactic space velocities U, V, W of our sample 
stars. These velocities are measured with respect to the Sun. The 
V component is measured in the direction of Galactic rotation, 
W is the component perpendicular to the plane of the Galaxy 
(positive towards the Galactic north pole), and U is measured 
towards the Galactic center. 

The transformation equations from radial velocity, parallax, 
and proper motion to Galactic space velocity we used are given 
by, e.g., Johnson & Soderblom (1987). Note, however, that due 
to correlations between astrometric parameters, as given in the 
Hipparcos catalog, the error estimates in this procedure need to 
be corrected (see the Hipparcos catalog or the Appendix B in 
AP04a). 

The typical error in the Hipparcos proper motion is about 
1 mas. Thus, the error in the heliocentric tangential velocities 
(V,) of our sample stars due to this error is very small. It is only 
about 0.05 km s _1 for a star at 10 pc and, in our sample, it reaches 
a maximum of 0.8 km s _1 for a distance of 150 pc. Therefore, the 
error in V, is dominated by the error in the parallax. For example, 
for the stars in our sample for which the error in the parallax is 
about 1%, the error in V, is about 0.5 km s _I while for the stars 
with a 10% error in the parallax the error in V, is about 10 km s _1 . 
Examining our sample results, we find that for thick-disk (thin- 
disk) stars, the mean error in V, is about 8.4 (1.6) km s _1 . Mean 
errors in U, V, W in our sample are 3.1, 4.0, and 2.5 km s _1 , re- 





Pi 


V 




cr v 


o~w 


i=l , thin disk 


0.90 


-12 


39 


20 


20 


i=2 , thick disk 


0.08 


-51 


63 


39 


39 


i=3 , halo 


0.02 


-199 


141 


106 


94 



Table 2. Parameterization adopted for the criterion of thin-disk, 
thick-disk, or halo membership. Units are km s _1 , except for p„ 
which is dimensionless. 



spectively. About 98 (91) % of our sample stars have uncertain- 
ties in V r smaller than 5(1) km s _1 which implies that the errors 
in V, dominate the uncertainties in the U, V, W velocities. For the 
nearest stars, the error in V r is normally comparable with the 
uncertainty introduced by the Hipparcos data, and only in a few 
extreme cases where the error in V r exceeds 10 km s _I , the errors 
in U, V, W are dominated by the uncertainty in V r . 

3.3. Kinematic criterion for subsampling 

Studies of stellar kinematics show that the distribution of U, V, W 
velocities may be reasonably well decomposed by assuming 
that there are three stellar populations involved, each with a 
Gaussian velocity distribution: thin disk, thick disk, and halo 
(e.g., Soubiran 1993, Soubiran et al. 2003). As a whole, the thick 
disk rotates slower around the center of the Galaxy compared to 
the thin disk and it has a wider velocity dispersion. The halo, 
on the other hand, has little or no overall rotation but the widest 
velocity dispersion. 

Typical values for the rotational lags with respect to the Sun 

are V 10 km s _1 for the thin disk, V 50 km s _1 for the 

thick disk, and V 200 km s _1 for the halo. There is still con- 
troversy regarding the exact values of both the rotational lags 
and the velocity dispersions but it is clear that the absolute value 
of the rotational lag and the velocity dispersions for the halo 
and thick disk are larger than those of the thin disk (compare, 
e.g., the value for the rotational lag of the thick disk by Nissen 

1995, 50 km s" 1 , with that given by Chiba & Beers 2000, 

~ -30 km s _1 , and the relatively large velocity dispersions for 
the thin disk according to Soubiran et al. 2003 with those given 
by AP04a)|]The mean thin/thick disk kinematics we adopted are 
given in Table [2] and correspond to those inferred by Soubiran 
et al. (2003). For the halo kinematics we adopted the values 
derived by Chiba & Beers (2000) for their 'pure' halo compo- 
nent. Note that the values in Table|2]are given in the heliocentric 
frame. 

Following the formulation adopted in previous studies (e.g., 
Mishenina et al. 2004, Bensby et al. 2005, R06), the probability 
that a star belongs to the component that has a mean V = V,- and 
a velocity dispersion cry = cr Vi , based only on the V velocity of 
the star, is proportional to exp[-(V - V,) 2 /(2o"y)]. Combining 
the three velocities, the probability of a star to belong to one of 
the components is 



(2tt) 3/2 c 



exp{-0.5 



U 2 (V 

2 2 



2 



°o-u,o- v ,o-w i 

where i = 1 for the thin disk, i = 2 for the thick disk, i 
the halo, and c, is a normalization constant given by 

Pi 



3 for 



Ci = 



zii Pi(Pi/cd 



7 However, it is possible that these differences are due to the use of 
samples with different characteristics. Note, for example, that the V ve- 
locities depend on spectral type (e.g., Dehnen & Binney 1998). 
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Fig. 2. Distribution of our sample stars in (U, V, W) Galactic space velocities. Crosses: Pi > 0.5 (thin-disk stars), circles: P 2 > 0.5 
(thick-disk stars), stars: P3 > 0.5 (halo stars). The dotted lines intersect at U = V - W = 0. 



In the last equation, the p, values are the relative number den- 
sities of thin-disk, thick-disk, and halo stars in the solar neigh- 
borhood. The equations are given in this form because the ratio 
Pi/cj can be computed directly from the observed data. 

The pi values are poorly known; values reported in the lit- 
erature for p2 (the relative density of thick-disk stars), for ex- 
ample, range between 0.02 and 0.15 (e.g., in order of increasing 
p 2 : Reid & Majewski 1993, Ojha 2001, Reyle & Robin 2001, 
Soubiran et al. 2003). Here we adopted p x = 0.90, p 2 = 0.08, 
P3 = 0.02. As revealed by star counts, the mean densities of 
thin-disk, thick-disk, and halo stars vary with distance from the 
Galactic plane (z). In a more general formulation, the pi's would 
need to be considered as functions of z. However, since our sam- 
ple stars are all nearby, the z-dependence does not have a notice- 
able effect and the use of constant local p, values is valid. 

Figs.|2]and[3]illustrate the kinematic properties of our sam- 
ple stars and the membership criterion adopted for subsampling. 
Table 4, available at the CDS, provides the basic data and kine- 
matic properties of our sample stars, including the probabilities 
described above. 

It is important to notice that the errors in the P, probabilities 
due to errors in the U, V, W velocities are only important for stars 
in the 'transition' regions of the Galactic velocity space. This is 
because even though there is an overlap in the kinematical distri- 
butions of the three groups, the peaks in V are still well defined 
(see Figs. [2] [3] and Table |2]i and a given error in the U, V, W set 
leads to a larger error in P, if P, ~ 0.5 compared to the P, ~ 1 .0 
case. By propagating the errors in the U, V, W velocities to the P, 
values of each star we find there are only 0, 3, and 1 thin-disk 
(with Pi > 0.7), thick-disk (P 2 > 0.7), and halo (P 3 > 0.7) 
stars, respectively, for which the uncertainty in P, is larger 
than 0.2 (the largest error found was 0.26). The corresponding 
mean values are < cr(P{) >= 0.005, < cr(P 2 ) >= 0.043, and 
< cr(P 3 ) >= 0.043. Therefore, if we set a P, > 0.7 constrain to 
determine the membership of the sample stars to the /-group, the 
chance of having a star with P, < 0.5 after considering errors in 
the U, V, W data is essentially zero. 



4. Abundance Analysis 

Both the iron and oxygen abundances were determined by 
matching the observed line equivalent widths with those pre- 
dicted by spectrum synthesis. A differential analysis was per- 
formed using the solar spectrum as a reference. 
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Fig. 3. Toomre diagram for our sample stars. Symbols are the 
same as in Fig. [2] Error bars are shown only for stars with un- 
certainties larger than 20 km s~' in Vf/ 2 + V 2 + W 2 . 



The classical 1D-LTE model atmospheres, without con- 
vective overshooting, by Kurucz (1993), interpolated to the 
(T e ff ,logg, [Fe/H]) set of each sample star, were adoptedQ 
Even though the solar model with convective overshooting re- 
produces better the solar spectrum, the derivation of stellar pa- 
rameters seems to be more robust with no-overshoot models 
(Castelli et al. 1997). Thus, we preferred the later type of Kurucz 
models for this work. 

The latest version of the spectrum synthesis code MOOG 
(Sneden 1973f| was used for the LTE calculations. Since the 
oxygen lines require a non-LTE treatment, the equations of sta- 
tistical equilibrium were solved with TLUSTY and non-LTE 



8 These are the models labeled 'nover' in the R. L. Kurucz's web- 
site: http : //kurucz . harvard . edu/gr ids . html We did not use the 
more recent and improved 'odfnew' models but we found that the 
changes in the abundances if we used them instead of the 'nover' mod- 
els are small (less than 0.02 dex in the case of iron) and in any case 
systematic so the relative abundances and abundance ratios are not af- 
fected. 

9 http : // verdi . as . utexas . edu/moog . html 
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spectrum synthesis was done with SYNSPEC (see Hubeny 1988, 
Hubeny & Lanz 1995)03 A detailed discussion of the non-LTE 
calculations is given separately in Sect. [5] 



4. 1 . Equivalent width measurements 

For almost all the spectral lines that we used, the Lorentzian 
wings could be neglected (or were avoided, since we did not use 
very strong lines affected by saturation) so the equivalent widths 
(EW) were measured by fitting Gaussian profiles to the spectral 
lines. The EW measurements were made using an automated al- 
gorithm that finds the central wavelength {A) of a line within a 
certain wavelength range and fits a Gaussian profile centered on 
this wavelength. Possible blends are taken into account in this 
procedure. 

The estimated error in our automatic EW measurements is 
3.5% (1.8 mA for a typical 50 mA line), based on a comparison 
of EW measurements of Fe i lines on the same star but using 
spectra from three of the data sources described in Sect. 12.11 
This comparison also showed no systematic differences in the 
EWs measured in spectra from different sources. 

The equivalent widths measured in this way for the three 
available solar spectra (see Sect. 12. It agree at the 1 % level, with 
the small differences in EW being random and not systematic. 
The adopted solar equivalent widths are a weighted mean of the 
EWs measured in each of the three spectra, with a higher weight 
given to the skylight spectrum from the S 4 N database. 



4.2. Line data and solar analysis 

We selected 119 Fei and 13 Fen lines to determine the [Fe/H] 
values of our sample stars. All these lines have transition proba- 
bilities measured in the laboratory and were chosen such that no 
evident blends with other features are present in the solar spec- 
trum. Our Fe i line selection also required the final line list not to 
show any correlation between the excitation potential (EP) and 
the reduced equivalent width (REW = EW/A) of the lines, which 
makes it ideal for purely spectroscopic studies, where these cor- 
relations lead to degenerate solutions for the stellar parameters 
that are estimated by making the line-by-line abundances inde- 
pendent of EP and REW. Note, however, that our procedure is 
not purely spectroscopic. 

We adopted laboratory gf- values for the Fe i lines from dif- 
ferent sources, all of which are described and compared to each 
other by Lambert et al. (1996). When available, the Oxford 
values (e.g., Blackwell et al. 1976) were adopted, otherwise a 
weighted mean with the typical uncertainties of the other sources 
was used. Transition probabilities for the Fe n lines are from the 
compilations of Allende Prieto et al. (2002) and Melendez et al. 
(2006). No astrophysical ^/-values were used. 

Oxygen abundances were determined only from the triplet 
O i lines at 777 nm. The transition probabilities adopted are those 
given in the NIST databaseQi.e., the average of the values cal- 
culated by Hibbert et al. (1991), Butler & Zeippen (1991), and 
Biemont & Zeippen (1992). In general, in order to analyze the 
[ O i] lines at 630 nm and 636 nm we would require higher resolu- 
tion and S/N. The exceptions are most of the S 4 N spectra, which 
are suitable for this type of analysis. However, we also avoided 
the 630 nm line due to the well-known blend with a Ni i line 
(e.g., Allende Prieto et al. 2001), which requires an accurate de- 
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Fig. 4. Our derived solar iron abundances from Fe i (circles) and 
Fen (squares) lines as a function of excitation potential and 
reduced equivalent width. Linear fits to the Fei data are also 
shown. 



termination of the Ni abundance prior to the oxygen abundance 
analysis. 

Van der Waals damping constants were adopted mainly from 
the Barklem et al. (2000) and Barklem & Aspelund-Johansson 
(2005) calculations. For the few Fei lines that are not included 
in the Barklem et al. (2000) work, the classical Unsold formula, 
with an enhancement factor of 2, was adopted. Standard ap- 
proximations for the radiative (e.g., Gray 1992) and Stark (e.g., 
Cowley 1971) broadening constants, as coded in MOOG, were 
also adopted. 

The line-list and atomic data adopted in this work are given 
in Table 5, available at the CDS. 



For the Sun, we obtained Ap e = 7.52 + 0.08 from the Fe i 
lines, and Ap e = 7.51 ± 0.07 from the Fe n lines[3 A microturbu- 
lent velocity of 0.90 km s _! was chosen to make the abundances 
from the Fei lines independent of REW. As shown in Fig. [4] no 
significant correlation with the EP of the lines is present either. 

Our mean solar iron abundance of Ap e = 7.5 ±0.1 is consis- 
tent, within error bars, with most other determinations using both 
classical 1D-LTE models (e.g., Holweger et al. 1991) and 3D 
hydrodynamical simulations (e.g., Shchukina & Trujillo Bueno 
2001 ; Asplund et al. 2004, 2005). 

An LTE solar oxygen abundance of Aq = 8.85 is ob- 
tained from the triplet. When non-LTE corrected, as in Sect. [5] 
this value reduces to Aq = 8.72, in good agreement with 
most recent determinations of the solar oxygen abundance (e.g., 
Allende Prieto et al. 2001; Asplund et al. 2004, 2005; Melendez 
2004; Shchukina et al. 2005). The agreement is perfect if we 
consider only the oxygen abundance derived from the triplet, 
corrected for non-LTE and using one-dimensional model atmo- 
spheres, but note that granulation effects reduce this abundance 
to 8.66 (see Asplund et al. 2004). 



10 http://nova.astro.umd.edu 

11 |http : //physics .nist . gov /PhysRefData/ASD/ 
lines_form.html; see also Wiese et al. (1996). 



12 In this paper, all error bars given for the abundances are lcr errors, 
where <r is the standard deviation of the line-by-line abundances. 
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4.3. Atmospheric parameters 
4.3.1. Effective temperature 

The infrared flux method (IRFM) color-temperature calibrations 
by Ramirez & Melendez (2005b) were used to estimate the effec- 
tive temperatures of our sample stars. They provide calibrations 
for 17 colors in different photometric systems. We used as many 
as possible of the 17 colors available for each star. On average, 
8 colors were used, independently of the population (thin-disk, 
thick-disk, or halo). 

The IRFM temperatures have been shown to be in good 
agreement with those derived from basic principles, i.e. with di- 
rect temperatures inferred from angular diameter and bolomet- 
ric flux measurements (e.g., Ramirez & Melendez 2004, 2005a). 
When compared to temperatures spectroscopically derived from 
the excitation balance of Fe i lines, however, a systematic differ- 
ence of about 100 K is observed for metal -rich main sequence 
stars (e.g., Ramirez & Melendez 2004, Santos et al. 2004, Yong 
et al. 2004, Heiter & Luck 2005). Interestingly, Santos et al. 
(2005) have found reasonably good agreement between spec- 
troscopic and direct temperatures for four of their sample stars 
while Casagrande et al. (2006) have shown that the zero point 
of the IRFM T e ff scale can be in good agreement with the spec- 
troscopic one depending on the flux calibration adopted. Masana 
et al. (2006) have also found a difference of 100 K with the IRFM 
T e g scale using their spectral energy distribution fit method. 
Besides this evidence for a systematic error in the zero point 
of the effective temperature scale, the use of several of the IRFM 
T e g vs. color calibrations by Ramirez & Melendez (2005b) guar- 
antees a very homogeneous, internally accurate, T e s scale for our 
sample stars. 
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Fig. 5. Difference between the [Fe/H] values obtained with the 
Fe ii and Fe i lines as a function of atmospheric parameters. Giant 
stars (logg < 3.8) are shown with triangles and dwarf stars 
(logg > 3.8) with squares. 



vs. Fe i difference. On the metal-poor side, Ramirez et al. (2006) 
have demonstrated that the spectroscopic \ogg value of at least 
one star is in disagreement with the more reliable surface gravity 
obtained from fitting of the wings of strong lines. 



4.3.2. Surface gravity 

The high accuracy of the Hipparcos trigonometric parallaxes al- 
lowed us to determine reliable surface gravities for our sample 
stars. The procedure is outlined in R03 and AP04a. 

An estimate of the stellar masses is needed to obtain the log g 
values. These were derived from a comparison of the location 
of the stars in the HR diagram with the theoretical isochrones 
by Bertelli et al. (1994). This approach also provides age esti- 
mates (Sect. [7). Errors in the mass determined in this way may 
be as high as 10%, considering not only errors in the fundamen- 
tal stellar parameters but also systematic errors in the theoretical 
isochrones. Nevertheless, this error affects only marginally the 
log g value whereas the error in the parallax affects it substan- 
tially (see, e.g., Eq. 2 in Allende Prieto et al. 1999). 

Given that the errors in the parallax of our sample stars 
are usually lower than 10%, the typical error in \ogg is about 
0.05 dex, which has only a marginal impact on the abundance 
determination. 

An alternative approach to the log g value is to tune it un- 
til the mean Fe abundances obtained separately from Fei and 
Fen lines agree, keeping all the other parameters constant. 
Allende Prieto et al. (1999) provide a comprehensive compar- 
ison of these spectroscopic surface gravities with those deter- 
mined from Hipparcos parallaxes. The spectroscopic procedure, 
however, is challenged by model uncertainties. Shchukina & 
Trujillo Bueno (2001), for example, have shown that overion- 
ization due to near-UV radiation may lead to a Fe n/Fe i ratio 
that differs from that obtained from the Saha equation in the so- 
lar photosphere. Similar results have been obtained by Thevenin 
& Idiail (1999) and Korn et al. (2003) but with a smaller Fen 



4.4. Abundances 

We determined absolute abundances from every spectral line of 
each star. Then, by comparing these results with the solar abun- 
dances, on a line-to-line basis, we inferred differential abun- 
dances for both iron and oxygen. 

The atmospheric parameters and iron abundances were de- 
termined iteratively starting with [Fe/H] = models. After three 
iterations, neither T e s nor log g showed further changes. At every 
iteration, a microturbulent velocity v, was determined by making 
the abundances from the Fe i lines independent of REW. 

AAA. The ionization balance problem 

The differences between the [Fe/H] values derived with Fe i and 
Fe ii lines are shown in Fig. [5] as a function of atmospheric 
parameters. It is clear that, on average, for T e $ > 5500 K, 
logg < 4.4, and [Fe/H] < 0.1, the Fen abundances are sys- 
tematically larger than the Fe i abundances by 0.06 dex. The dis- 
agreement worsens for cooler, higher surface gravity, and more 
metal-rich dwarf stars. Cool giants with 4700 K< T eff < 5000 K 
also show systematically large Fe n abundances but at a roughly 
constant level of about 0.06 dex. The solution to this problem is 
at present unknown, although it could be explained by system- 
atic errors in the stellar parameters, the inadequacy of classical 
abundance analysis that use static one-dimensional model atmo- 
spheres and LTE line formation calculations, or a combination 
of them. 

For most of our sample stars, systematic offsets in the T e g 
and log g scales would make the Fe i and Fe n abundances agree. 
For a solar- like star, for example, an increase of 100 K in T e g 
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leads to an increase of 0.07 dex in the Fe i abundance and a de- 
crease of 0.03 dex in the Fen abundance. If the surface grav- 
ity is then increased by 0. 1 dex the Fe i abundance decreases by 
0.01 dex while the Fen abundance increases by 0.04 dex. Thus, 
after an increase of 100 K in r e jf and 0.1 dex in logg, the Fei 
abundance increases by 0.06 dex while leaving the Fe n abun- 
dance nearly unchanged. Excluding the cool metal-rich dwarfs, 
this would 'solve' the ionization balance problem shown in 
Fig. 13 

Interestingly, a systematic difference of 100 K between the 
IRFM and the excitation balance temperature scales has been re- 
ported in the literature (Sect. l4~3~TT >. Also, on average, the ioniza- 
tion equilibrium gravities seem to be larger than those obtained 
using Hipparcos parallaxes by roughly 0. 1 dex (Allende Prieto 
et al. 1999). Note that if the effective temperatures are 100 K hot- 
ter, the trigonometric logg values would increase by 0.07 dex 
(cf. eq. 2 in Allende Prieto et al. 1999) so there is no need to 
question the accuracy of the Hipparcos parallaxes in any case. 

Alternatively, if the stellar parameters we adopted are cor- 
rect, we should question the simplifications made in the abun- 
dance analysis, in particular the use of static plane parallel model 
atmospheres and LTE line formation. A full discussion of the 
effects of surface inhomogeneities and non-LTE effects on the 
Fe line formation is beyond the scope of this paper (see, e.g., 
Thevenin & Idiart 1999, Asplund et al. 2000, Shchukina & 
Trujillo Bueno 2001, Korn et al. 2003, Shchukina et al. 2005). 
Nevertheless, evidence suggests that non-LTE effects may be re- 
sponsible for the Fe i vs. Fe n discrepancy. 

Ionization by UV photons decreases the number density of 
neutral Fe atoms and increases that of Fe n ions. In warm stars 
(Jeff > 5500 K), most of the Fe is in the ionized form so the 
main effect on the abundances is the decrease of the Fe i density. 
Thus, in a statistical equilibrium calculation where overioniza- 
tion by UV photons is the main non-LTE effect, the Fei lines 
are weaker than those predicted in LTE. In that case, to match 
the observations, the Fe i abundance has to be increased com- 
pared to the LTE case, thus reducing the difference with the Fe n 
abundances, which do not change significantly due to this ef- 
fect. Note, however, that this is expected to be less important in 
metal-rich cool stars, which contradicts the fact that the Fe i vs. 
Fe ii discrepancy worsens for our K-dwarfs of solar and super- 
solar metallicities (Fig. [5]). Other model uncertainties, for exam- 
ple granulation or activity-related effects, may account for these 
large differences. Also, due to line blanketing, the overionization 
effect described above may result in a stellar atmosphere with a 
temperature structure different from that obtained in LTE, which 
may be also responsible for the differences seen at high [Fe/H], 
as suggested by Feltzing & Gustafsson (1998). 

4.4.2. Adopted iron abundances 

Due to the ionization balance problem described above, instead 
of averaging the abundances from Fe i and Fe n lines, we would 
prefer to use the mean abundances of iron from the Fe n lines 
only as our [Fe/H] scale due to their relatively low sensitivity 
to the atmospheric parameters adopted for most of our sample 
stars, i.e. those hotter than 5400 K, and small predicted non-LTE 
effects (e.g., Thevenin & Idiart 1999, Gehren et al. 2001). 

Unfortunately, errors in the atomic data and the lower num- 
ber of Fe ii lines employed in this work compared to the number 
of Fei lines make the Fei abundances (most likely incorrect) 
internally more accurate. In fact, the mean line-to-line scatter of 
the Fe i abundances is 0.06 dex while that of the Fe n abundances 
is 0.10 dex. Thus, instead of discarding the Fe i abundances, we 
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panel: T e g vs. [Fe/H], as derived from the Fen lines. Lower 
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Crosses: P\ > 0.5 (thin-disk stars), circles: P2 > 0.5 (thick-disk 
stars), stars: P3 > 0.5 (halo stars). 



prefer to embrace the scale of the Fe n abundances by increas- 
ing the Fe 1 abundances by 0.06 dex. The adopted [Fe/H] value 
is then a weighted mean of the modified Fe 1 abundance and the 
Fe n abundance, where the weights are the individual values of 
the line-to-line scatter. Since the 0.06 dex difference is not con- 
stant with stellar parameters, we applied this correction only to 
stars with T e ff > 5400 K, logg > 3.7 (extrapolation to higher 
metallicities was allowed so the reader should be careful when 
using our [Fe/H] values for stars with [Fe/H] > 0.1). 

The adopted atmospheric parameters and derived iron and 
oxygen abundances are given in Table 6, available at the CDS. 
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Figure [6] illustrates the overall properties of our sample in the 
HR diagram, as well as abundance vs. r e ff diagrams. 

5. Non-LTE corrections for the O i triplet 

The O i triplet lines are known to be formed under conditions 
far from the LTE approximation, both in giant (Eriksson & Toft 
1979) and main sequence stars (Kiselman 1993, 2001). Our 
non-LTE corrections were calculated using fixed temperature 
and electron density structures from Kurucz LTE model atmo- 
spheres. Only the level populations of the hydrogen and oxy- 
gen atoms were allowed to depart from LTE by solving the rate 
equations and consistently recalculating the radiation field with 
TLUSTY. 

We used an oxygen model atom with 54 levels plus the con- 
tinuum and 242 transitions. Details on the atomic data used in the 
non-LTE calculations are given by Allende Prieto et al. (2003a). 
In particular, note that the model atom does not consider the fine 
structure of the triplet. Allende Prieto et al. (2003b) used this 
model atom to show that the lower level of the triplet is slightly 
overpopulated with respect to the LTE case while the upper level 
is noticeably underpopulated in the solar case. They also showed 
that the levels corresponding to the [Oi] 630 nm line do not 
depart substantially from LTE. These results are fairly consis- 
tent with those obtained by other authors (e.g., Kiselman 1993, 
Takeda 1994, Gratton et al. 1999). 

Once our calculations were complete, we discovered an error 
in the collisional ionization ratesQ Nevertheless, the difference 
between the predicted equivalent widths (EW) of the triplet lines 
computed using the old and new model atoms is small. In gen- 
eral, the predicted EWs decrease by about 3% and only a small 
dependence with temperature and metallicity is observed. The 
predicted EWs decrease by 3.2% for a cool (T e g = 5000 K) star 
and by 2.4% for a hot star (T e ff = 6400 K) while they decrease 
by 2.8% in the solar case and by 3.8% for a [Fe/H] = -0.8 star 
of solar temperature. In all cases, the effect on the oxygen abun- 
dance (Ao) derived by comparing the predicted EWs with the 
observed ones was an increase of 0.02 dex (rounded to the sec- 
ond decimal point). Therefore, the [O/H] ratios are unaffected 
given that the same abundance shift needs to be applied to our 
reference star (the Sun). 

The non-LTE corrections we calculated are shown in Fig. [7] 
They are larger for warmer and low surface gravity stars in which 
the radiation field is stronger and the densities lower thus in- 
creasing the importance of radiative transitions and decreasing 
that of collisional excitation. More elegant descriptions of the 
non-LTE effects involved in the formation of the triplet can be 
found in Eriksson & Toft (1979) and Kiselman (1993). 

An important but still uncertain ingredient for the non-LTE 
corrections are inelastic collisions with neutral H atoms. If in- 
cluded, the non-LTE corrections will be reduced. Allende Prieto 
et al. (2004b) have shown that H collisions, although hav- 
ing a marginal and usually neglected role at solar metallicity 
(e.g., Nissen et al. 2002), help to reproduce better the center- 
to-limb variation of the triplet line profiles when included using 
Drawin's formula (Drawin 1968) multiplied by an empirical fac- 
tor (S h), as suggested by Steenbock & Holweger (1984). Note, 
however, that this effect becomes more important at low metal- 
licities (Ramirez et al. 2006). 

The key effect of introducing H collisions is the decrease 
of thermalization with depth due to the decrease of neutral H 

13 A corrected model atom has been made available at 
http : / /hebe . as . utexas . edu/ at 



density, which leads to smaller non-LTE corrections for the 
strongest line of the triplet compared to those for the weakest 
line (see Allende Prieto et al. 2004b and Ramirez et al. 2006). 

Our non-LTE calculations do not include the effect of H col- 
lisions. In this way, by analyzing our results we find that in 
order to obtain consistent non-LTE abundances from the three 
lines of the triplet, the non-LTE corrections of the 7771.9 A and 
7774.2 A lines have to be decreased by 0.036 dex and 0.018 dex, 
respectively, to match the abundance from the 7775.4 A line 
(nearly independent of stellar parameters). These mean empir- 
ical corrections were applied to force an agreement between the 
abundances derived from each line of the triplet, which in turn 
reduced the scatter in the oxygen abundance patterns but may 
be still giving systematically low abundances, with, perhaps, a 
trend with temperature and metallicity. Note that, for example, 
Allende Prieto et al. (2004b) find an increase in the solar O abun- 
dance of 0.06 dex when going from a non-LTE calculation with- 
out H collisions (S h = 0) to one with them (Sh = 1) while 
Ramirez et al. (2006) find an increase of 0.25 dex for a star with 
[Fe/H] = -1.7 when going from S H = to S H = 100 

For the Sun, we obtained a mean non-LTE correction of 
0.13 dex, which leads to a non-LTE oxygen abundance of Ao = 
8.72 + 0.02, in agreement with other recent solar oxygen abun- 
dance determinations, as mentioned in Sect. 14.21 

As in the case of the iron abundances, for the sample stars, 
oxygen abundances were calculated differentially on a line to 
line basis with respect to the Sun, and are given as [O/H] in 
Table 6, available at the CDS. 

5.1. A grid of non-LTE corrections 

The non-LTE abundances given in Table 6 were computed 
'directly' from curves of growth corresponding to the stel- 
lar parameters of each particular star. For practical purposes, 
we also computed a grid of non-LTE corrections in the 
r e jf,logg, [Fe/H],Ao,L space, which allows to obtain non-LTE 
corrections, given the atmospheric parameters and LTE oxygen 
abundances from each line of the triplet, for any star with param- 
eters within those covered by our sample. The grid and an IDL 
routine to interpolate within this grid are available from the au- 
thors upon request. The non-LTE abundances obtained directly 
and those inferred using the grid agree at the 0.02 dex level with- 
out any trends with stellar parameters. 

6. Oxygen Abundance Patterns 

Here we describe the oxygen vs. iron abundance trends we de- 
rived for halo, thin- and thick-disk stars using the preceding re- 
sults. Interpretations of the results are given in Sect. [8] 

Fig. [8] shows the trends we obtained before and after apply- 
ing the non-LTE corrections to all our sample stars. It is clear 
that the non-LTE corrections reduce the scatter and make the 
thin/thick disk differences obvious compared to the LTE case. 
In the LTE case, the [O/Fe] vs. [Fe/H] relation shows a steeper 
slope compared to the non-LTE case. Also, in the LTE case, no 
clear differences are seen between thin- and thick-disk stars, ex- 
cept in the [Fe/O] vs. [O/H] plane, where an apparent continuous 
increase in the [Fe/O] ratios from thick- to thin-disk stars is seen 
at all [O/H]. The fact that the thin/thick disk separation becomes 
obvious after applying the non-LTE corrections suggests that our 
non-LTE calculations are reliable. 

14 In each case, the S h + values were chosen to reproduce accu- 
rately the triplet line profiles. 
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Fig. 7. Difference between the oxygen abundance derived using LTE (Ao,l) and non-LTE (Ao,n) level populations, i.e. non-LTE 
corrections, for the middle line of the 777 nm O i triplet as a function of atmospheric parameters and equivalent width. 
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Fig. 9. As in the right side panels of Fig.[8]but restricted to stars with T e g > 5400 K, logg > 3.8, and membership probabilities 
larger than 70%. The [Fe/H] values here are a mean of the Fe n and modified Fe i abundances (Sect. 14.4.21 ). Squares and triangles 
represent stars whose kinematic probabilities of being either thin-disk, thick-disk, or halo stars are lower than 70%. Triangles are 
stars that, kinematically, may belong to either the thick disk or halo, while squares represent stars that may belong to the thin disk 
or the thick disk. The dotted, solid, and dashed lines are linear fits to the thin-disk, thick-disk, and halo samples, respectively (see 
Table [3), when the kinematic probabilities are larger than 70%. For the thick-disk fit, we excluded the stars that having thick-disk 
kinematics follow the thin-disk abundance pattern. 



To reduce the impact of errors in the atmospheric parameters 
and model uncertainties, we also show the oxygen abundance 
trends for stars with > 5400 K and logg > 3.8 (i.e. dwarfs 
only) in Fig. [9] In this way, besides restricting the sample to stars 
with similar parameters, we also avoid the coolest dwarfs that 
suffer more strongly from the ionization balance problem dis- 
cussed in Sect. 14.4. 71 In fact, the [Fe/H] values in Fig. |9]corre- 
spond to those derived in Sect. 14.4.21 where the Fe i abundances 
are empirically put into the more reliable Fe n abundance scale. 
Furthermore, in Fig.|9]we discriminate between stars with high 
and intermediate probabilities of being thin-disk, thick-disk, or 
halo stars by adopting the strong constrain of P, > 0.7, which 
rules out any possibility for the errors in the U, V, W velocities of 
a given star to allow Pj to be lower than 0.5 (see Sect. |3.31 >. 

6.1. Halo, thin-disk, and thick-disk stars 

We computed linear fits to the data shown in Fig. [9] separating 
the stars into halo, thin-disk, and thick-disk members using the 
strong constraint of 70% in the kinematic membership probabil- 
ities. Also, the thick-disk fit does not include the kinematically- 
selected thick-disk members that follow the thin-disk abundance 
pattern. These are discussed separately in Sect. I6.3I The results 
of the fits are given in Table [3] For these calculations, the ran- 
dom errors in [Fe/H] and [O/Fe] for each star were taken into 
account. 

Although there are a few very metal-poor halo stars within 
150 pc from the Sun, we restricted our sample of halo stars to 



[Fe/H] > -1.5 because our procedure, specifically the auto- 
matic measurement of equivalent widths (Sect. 14. It , works best 
at these metallicities. Therefore, our halo star sample is not com- 
plete, and any attempt to trace a halo abundance trend using only 
the metal-rich end is dubious. Nevertheless, there is an overlap in 
the [Fe/H] distributions of halo and thick-disk stars from about 
-1.2 to -0.4 which can be analyzed using the results from this 
study. In this range the mean [O/Fe] ratios of thick-disk and halo 
stars are similar at about 0.45. By inspecting the linear fits in 
Table [3] however, we notice that only the thick disk trend shows 
a significant negative slope. Halo stars seem to have, statistically, 
constant [O/Fe] in this metallicity range while the [O/Fe] ratio 
of thick disk stars decreases by almost 0. 1 dex. Note also that 
the star-to-star scatter for the halo star sample is slightly larger 
than that of the thick disk. 

A quick inspection of Figs. [8] and [9] reveals that the ma- 
jority of thick-disk stars have larger [O/Fe] and [O/H] ratios 
compared to thin-disk stars for [Fe/H] < -0.2. Ignoring the 
thick-disk stars that appear to follow the thin-disk abundance 
trend, we conclude from Figs. [8] and [9] that there are two dif- 
ferent and mostly separated abundance patterns for thin- and 
thick-disk stars. Compared to thin-disk stars, thick-disk mem- 
bers have larger [O/Fe] and [O/H], and smaller [Fe/O] in the 
ranges -0.9 < [Fe/H] < -0.2, -0.5 < [O/H] < 0.3, respec- 
tively. The differences are roughly constant at about 0.2 dex in 
[O/H] and 0.3 dex in [Fe/O], but they depend on [Fe/H] in the 
case of [O/Fe] because the thin disk [O/Fe] ratio has a steeper 
slope (Table [3}. Note also that while the thick disk trend is well 
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a 


b 




n 


0" 


Thin disk 
Thick disk 
Halo 


0.096 
0.370 
0.388 


± 0.004 
± 0.027 
± 0.049 


-0.327 ± 
-0.121 ± 
-0.048 ± 


0.016 
0.043 
0.071 


238 
43 
12 


0.056 
0.065 
0.072 



Table 3. Linear fits to the data shown in Fig. [9] Given here are the 
coefficients of the fit: [O/Fe] = a + b[Fe/H], the number of stars 
used («) and the random scatter around the mean fit (<x). The 
linear fits and the cr computation take into account the individual 
errors in [Fe/H] and [O/Fe]. 



fitted by a line, the thin disk trend seems to require a higher order 
polynomial to satisfactorily fit the data. 

We estimate the random error due to mean errors of 50 K in 
r e ff and 0.05 dex in \ogg to be about 0.03 dex for [Fe/H], the 
line-to-line scatter is 0.06 for [Fe/H] and 0.03 for [O/H] after 
correcting for non-LTE, so the random error in [O/Fe] should be 
about 0.07 dex. We expect this random error to be slightly larger 
for the halo sample due to uncertainties in reddening, parallax, 
and lower S IN spectra. It is also reasonable to conclude that the 
scatter is smaller for the thin-disk sample, for which better spec- 
tra are available, as well as accurate stellar parameters. Thus, 
we find that there is little or no room for cosmic scatter in the 
thin-disk and the thick-disk abundances, in agreement with the 
results of R03 and R06. Although we can be tempted to draw the 
same conclusion for the halo sample, we should be cautious and 
warn the reader that, before that, more metal-poor stars need to 
be added to the sample. Note that, for example, Melendez et al. 
(2006) find that the observed scatter in the [O/Fe] ratios of halo 
stars in the -3 < [Fe/H] < -1 range allows for a small intrinsic 
scatter. 

6.2. An intermediate population? 

Our sample is biased towards stars with high kinematic probabil- 
ities of being either thin- or thick-disk stars because it is domi- 
nated by nearby stars (hence the large number of thin-disk mem- 
bers) and a sample selection was specifically made to collect a 
significant number of stars with kinematics similar to those of 
the thick disk. Thus, only a few stars (18) in our sample have 
low kinematic probabilities of being either thin- or thick-disk 
stars and only these should be looked at to check if an interme- 
diate population with properties between those of the two disks 
exists. These 18 stars are shown in Fig.fTOh. 

A quick inspection of Fig.fTOh (see also Fig.[9j» already shows 
that stars with a kinematics intermediate between that of the 
thin and thick disk seem to clearly follow either the thin disk 
or the thick disk abundance pattern and not an intermediate one. 
To show this quantitatively, we simulated the scatter around the 
mean trends for the stars in our sample that have kinematic prob- 
abilities of being thin- or thick-disk stars below 70% assuming 
that: (I) they belong to a population with an abundance pattern 
between those of the thin and thick disk (Fig.fTOb), and (II) they 
belong to either the thin disk or the thick disk (Fig. fTOb). The 
mean thin/thick disk trends and scatter adopted here are those 
given in Table[3] For the intermediate population we assumed an 
scatter of 0.06 dex, which is the average of the scatter observed 
in thin- and thick-disk stars. 

Clearly, case (I) provides a poor representation of the data. 
Even if we assume that the intermediate population has a larger 
scatter around the mean abundance trend, it would be hard to 
explain the dip in the scatter distribution shown in Fig.fTOb. On 
the other hand, case (II) satisfactorily reproduces the dip and 
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Fig. 10. (a) Stars in our sample with low kinematic probability 
(< 70%) of being either thin- or thick-disk stars in the [O/Fe] vs. 
[Fe/H] plane. The solid lines represent the mean trends of thin- 
and thick-disk stars with high kinematic probabilities in our sam- 
ple, as given in Table [3] The dotted line represents the interme- 
diate abundance trend of a ficticious population with chemical 
properties between those of the thin and thick disk, (b) Case I: 
simulated stars (18, the same number of stars as in (a)) assum- 
ing that they all belong to the intermediate population, (c) Case 
II: simulated stars assuming that half (9) belong to the thin disk 
and half (9) to the thick disk, (d) Scatter distribution of the ob- 
served (panel (a), thick solid line histogram) and simulated [Case 
I (panel (b), dashed line), Case II (panel (c), thin solid line)] data. 
Gaussian fits are used for the simulated data. 



the peaks in Fig. [TUb even though the statistics is poor (only 
18 points). A Kolmogorov-Smirnov test (see, e.g., Press et al. 
1992) shows that the probability that the data shown in Fig.fTOb 
belong to the double peaked thin/thick disk distribution (the thin 
solid line in Fig.fTOb) is 96% while the significance level for the 
intermediate population case (the dashed line in Fig.fTOb) is less 
than 1%. Only if we assume that the scatter around the mean 
trend of the intermediate population is 0.14 dex, this number 
increases to a maximum of 72%. In such case, we would need a 
cosmic scatter of about 0.12 dex in the intermediate population. 
Given that no cosmic scatter is seen in the thin or thick disks 
(Sect. 16. IT ), this wide distribution seems unrealistic. This analysis 
discards the existence of a population of stars with abundances 
intermediate between those of the thin and thick disk, at least in 
the solar neighborhood. 
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6.3. The 'knee' and Reddy et al. TKTA stars 

In a series of papers, T. Bensby and collaborators have claimed 
that the abundance trends of the thick disk show the contribution 
of Type la supernovae (Feltzing et al. 2003; Bensby et al. 2003, 
2004, 2005). In the case of oxygen, Bensby et al. (2004, hereafter 
B04) found that from an almost constant [O/Fe] = 0.35 at low 
metallicities, thick disk stars with [Fe/H] > -0.3 show [O/Fe] 
ratios that decrease rapidly towards [O/Fe] = at [Fe/H] = 0, 
i.e. they find a 'knee' in the [O/Fe] vs. [Fe/H] relation of thick- 
disk stars. 

In Fig. [9] there is a significant number of stars with high 
kinematic probability of being thick-disk stars having [Fe/H] > 
-0.4. Some of these were recently observed by us and are thus 
not included in the Reddy et al. papers. We targeted this sample 
aiming at a better definition of the knee and/or determine where 
the thick disk [Fe/H] distribution ends. However, we found that 
almost all of these stars follow the thin disk abundance pattern. 
In fact, in Fig.[8]we also see some thick-disk stars at even lower 
metallicities that follow the thin disk abundance pattern. This 
suggests that the kinematics alone does not always correctly sep- 
arates thin- from thick-disk members. The possibility of errors in 
the abundance determination for these outliers is very unlikely. 
Therefore, following R06 designation (see their Sect. 6.1), these 
are TKTA stars (stars with Thick-disk Kinematics but Thin-disk 
Abundances). 

The TKTA stars are probably either old thin-disk stars that 
have been heated by secular processes (more than the average) 
or members of a moving group that can be ultimately associ- 
ated with the thin disk. As shown by R06, the TKTA stars are, 
on average, older than a typical thin-disk member but apparently 
younger than a typical thick-disk star and their distribution in the 
U, V, W galactic velocity space shows a preference towards pos- 
itive U velocities even though thick-disk members cover both 
positive and negative U velocities, more or less uniformly, at 
least in the ±150 km s _1 range. All the stars we classify as TKTA 
are older than 4 Gyr (Fig. [T3l but their distribution in U, V, W 
is more random compared to the R06 case which implies that 
only some TKTA stars (most of them in the R06 sample) may be 
members of a moving group. Note also that the Galactic bar may 
have affected the orbits of some low eccentricity stars in such a 
way that they now resemble thick-disk members, as suggested by 
Soubiran & Girard (2005). This may also explain the existence 
of stars with low kinematic probabilities of being thin- or thick- 
disk members, which, instead of being members of an interme- 
diate population, show clearly thin-disk abundances (Sect. l6T2l . 

Two thick-disk stars in Fig. [9] with [Fe/H] slightly below 
solar do show clearly higher [O/Fe] ratios compared to the 
mean thin disk (HIP 99224 and HIP 90542). They both have 
r e ff st 5500 K, closer to our lower limit of 5400 K, a temper- 
ature at which the triplet lines begin to give systematically large 
abundances for young cool dwarf stars (Morel & Micela 2004, 
Schuler et al. 2006), most likely due to the use of classical model 
atmospheres that do not include activity related and/or granu- 
lation effects. Note also that this is roughly the temperature at 
which the ionization balance problems described in Sect. 14.4.11 
begin to become important. If we restrict the sample to stars 
with r e ff > 5600 K, the 'oxygen-enhanced' thick disk ends at 
[Fe/H] = -0.3, again without obvious evidence for a knee (see 
Fig. ITTb . Therefore, our sample, expanded to better define the 
metal-rich end of the thick disk, shows no signature of a knee in 
the thick disk [O/Fe] vs. [Fe/H] relation. Given that others have 
proposed that metal-rich thick-disk stars exist and define a knee 
(in particular Bensby et al. 2004), our solid conclusion from this 
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Fig. 12. (a) LTE and non-LTE Abundance trends for the thin- 
(open circles) and thick-disk (filled circles) stars in Bensby et al. 
(2004) using their atmospheric parameters but our atomic line 
data and non-LTE corrections, (b). As in (a) but using our J e ff 
and \ogg scales. The dotted and solid lines are the mean trends 
found with our sample (TableO. 



is that not all metal-rich stars with thick disk kinematics show 
that behavior. In fact, the vast majority of them do not follow the 
knee, if there is any at all. 

6.4. Re-analysis of the Bensby et al. data 

The equivalent widths for the lines used in B04 were kindly sup- 
plied by T. Bensby. Using these data and the atmospheric pa- 
rameters from their paper, we obtained LTE iron and oxygen 
abundances, and corrected for non-LTE effects using our grid 
(Sect. [5TTT >. The results are shown in Fig.[T2h. 

Our LTE oxygen abundance trends, when using the B04 pa- 
rameters, are essentially the same as those obtained by B04 (see 
their Fig. 9c), as expected. The small differences are only due to 
the use of slightly different log gf values. When corrected us- 
ing our non-LTE calculations, no clear evidence of a knee is 
observed, in contrast with the B04 non-LTE result, which was 
obtained using non-LTE corrections from Gratton et al. (1999). 
Comparing their Fig. 9d with our Fig [T2"a we notice that it 
is essentially one star (HIP 109450: [Fe/H] ~ -0.15 in B04, 

[Fe/H] 0.10 in our work) that 'defines' the knee in the B04 

work (using the triplet, their knee obtained using the results from 
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Fig. 11. As in Fig. [9]but restricted to thin- and thick-disk stars with r e (j > 5600 K and membership probabilities larger than 70%. 
The dotted and solid lines correspond to the mean trends given in Table [3] using a larger sample, the thick dashed line shows the 
predictions of our simple GCE model for the thick disk (the details of the model are given in the Appendix [B] and the discussion 
in Sect. [8j, and the thick dotted line is the model by Nomoto et al. (2006) shifted upwards by 0. 1 dex to match our results at low 
metallicity. 



the forbidden 630 nm [Oi] line is defined by more stars). This 
star has a near solar temperature but it is slightly more evolved 
(log g ~ 4.2). Thus, its non-LTE correction should be larger than 
that of the Sun by about 0.1 dex (Fig. [7]) leading to a non-LTE 
[O/Fe] ratio smaller by roughly the same amount when com- 
pared to the LTE case. However, a visual inspection of Fig. 9d of 
B04 suggests that their LTE and non-LTE [O/Fe] ratios for this 
star are very similar. Therefore, the non-LTE corrections adopted 
by B04 may not be correctly taking into account the effects of the 
surface gravity. Our Fig.fTZh shows that the [O/Fe] ratio of this 
star has indeed decreased by 0. 1 dex when going from LTE to 
non-LTE. As a result, in our FigQUi, HIP 109450 is clearly fol- 
lowing the thin disk abundance pattern, contrary to the B04 case, 
in which the star is in between the thin- and thick-disk trends, 
thus creating an apparent knee. 

In Fig [T2b we show the abundance trends obtained with the 
B04 data but using effective temperatures, metallicities, and sur- 
face gravities computed as in our work. The systematic differ- 
ences in the abundances, compared to B04, are mainly due to the 
use of different temperature scales. When non-LTE corrected, 
the abundance trends derived from the B04 data are in excellent 
agreement with those found in our large sample (TableO, shown 
in Fig. [12b with solid lines. This shows that the differences with 
the B04 work for the triplet lines, in particular the discrepancy 
regarding the existence of the knee, are not a consequence of 
sample selection or data quality. 

Note that Bensby et al. (2005) add more stars to the B04 
sample, thus strengthening their O vs. Fe trends (including some 
'knee stars') obtained using the 630 nm [O i] line. Unfortunately, 
the spectral region where the triplet is was not included in their 
latest observations (T. Bensby, private communication) and thus 
no further comparisons are possible. 

7. Stellar ages 

To complement the observational evidence, we computed the 
ages of our sample stars using our derived stellar parameters and 
the Bertelli et al. (1994) isochrones, as explained in R03 and 
AP04a. These ages may be subject to large systematic uncer- 
tainties, as explained in the Appendix lAl but in order to main- 
tain consistency with the other stellar parameters derived from 
the same set of isochrones (namely the logg values), we did not 
correct for those. 



In Fig. [13] we plot the stellar ages vs. abundances for the 
stars in Fig. [9] Only the stars for which precise ages could be 
determined (<x(log (Age/yr)) < 0.3 dex, where cr is the internal 
uncertainty) are shownCjIt is clear that our thick-disk stars with 
[O/Fe] > 0.25 are older than our thin-disk stars and that the 
former seem to cover a short age range (10-14 Gyr) while thin- 
disk stars as young as 2.5 Gyr or less and as old as 10 Gyr are 
found. These ages are in reasonable agreement with those de- 
rived by Fuhrmann (1998), who suggests that all thin disk stars 
are younger than 9 Gyr while all thick disk stars are older than 
10 Gyr, and Gratton et al. (2000, their Figs. 5d and 6d). An in- 
teresting conclusion may be drawn from these results: the thick 
disk was formed before the thin disk, or at least the bulk of it. 
Whether this implies a hiatus of star formation, and therefore an 
age gap in Fig. [13] between the formation of the thick and thin 
disks (as suggested by Fuhrmann 1998 and Gratton et al. 2000 
based on the abundance information) depends on the accuracy 
of the ages. An age gap might be present, but the uncertainties 
in our derived ages make difficult to detect it. We remark that 
the majority of the young stars (Age < 8 Gyr) identified, kine- 
matically, as thick-disk members in Fig. [13] follow the thin-disk 
abundance pattern. 



8. Discussion 

Here we discuss the implications of our results for thin- and 
thick-disk stars. Our halo star sample is by far incomplete so 
we will not discuss the relevance of our results for halo studies. 

Oxygen is believed to be produced mainly in massive stars 
and put into the interstellar medium (ISM) by Type II supernovae 
(SN II) explosions, while iron is produced in both SN II and 
Type la supernovae (SN la) events. The time scale for SN II ex- 
plosions to occur is about 3 x 10 7 yr, which roughly corresponds 
to the lifetime of a star with a mass of 8M . For a SN la event, 
however, one must wait for a low or intermediate mass star to 



15 The lower and upper limits in the age of a given star were calculated 
from its age probability distribution (APD). The lower (upper) limit 
corresponds to the value at which the integrated APD, starting from 
the peak, reaches 0.955/2, i.e., the equivalent of the area covered by 
a Gaussian APD from the peak to the lower (upper) 2-cr limit. This 
procedure results in more realistic error bars given the asymmetry of 
the APDs. 
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Fig. 13. Age vs. abundance for the stars shown in Fig. [9] for 
which reasonably accurate ages can be derived. The dashed line 
corresponds to the age-metallicity relation adopted in our simple 
GCE model (AppendixlBl with the solid line showing the first 4 
Gyr of that relation. 



evolve into a white dwarf and accrete material from a compan- 
ion. Since lifetimes of these stars are substantially longer, the 
time scale for the first SN la events to take place is also longer, 
of the order of ~ 10 9 yr, although it should be noted that the ex- 
act value for the time scale of SN la events is very uncertain (see, 
e.g., Greggio 2005 for a recent discussion on this topic). In prin- 
ciple, stellar oxygen and iron abundances, as well as their ratios, 
contain information about the chemical enrichment of the ISM. 
This simple picture is complemented by additional ingredients 
such as the supernova yields that determine the net amount of 
the different elements ejected in the explosions, and the star for- 
mation rate and initial mass function (e.g., Wheeler et al. 1989, 
McWilliam 1997). 

Intermediate mass stars (those with masses between about 
1 and 10 solar masses) also contribute a significant fraction of 
metals to the ISM (e.g., Iben & Renzini 1983). During their evo- 
lution, in the so-called Asymptotic Giant Branch (AGB) phase, 
these stars undergo a series of thermal pulses that cause heavy 
element production (particularly C, N, and ^-process elements) 
in the helium burning shells, which are later ejected with the en- 
velope as winds and superwinds. 

Note that the stellar yields may strongly depend on the ini- 
tial metallicity and the minimum mass for black hole formation 
(e.g., Schild & Maeder 1985, Maeder 1992). In the former case, 
metallicity affects stellar evolution through changes in the opac- 
ities, which affect both the internal structure and mass loss rates. 
Oxygen production seems to be favored at low metallicities (e.g., 
Kobayashi et al. 2006). Regarding black hole formation, an un- 
certainty in the minimum mass required to form a black hole 
(Mbh) leads to an uncertainty in the net heavy element yield. 
If Mbh 1S large* then more stars (those with lower masses) will 
effectively eject their outer layers rich in heavy elements. On 
the other hand, if Mbh 1S small, more stars will become black 
holes and swallow part of the heavy elements during the col- 
lapse. Thus, the heavy element yields decrease with a decreas- 
ing Mbh- Note also that the net yield depends on both the initial 
metallicity and minimum mass for black hole formation, given 
that some of the synthesized material is ejected by winds. While 
this is true for SN II events that lead to black hole formation, in 



general the amount of metals returned to the ISM depends on the 
mass of the compact remnant, which may be also a neutron star. 

The abundance trends we derive for the thin and thick disks 
(Figs. [8] [9] and [TT|i suggest that the gas from which thick-disk 
stars were formed was rich in material synthesized in massive 
stars, given the larger [O/H] and particularly [O/Fe], essentially 
at all [Fe/H]. The thick disk [O/Fe] ratio, which is, in a first 
approximation, roughly constant at about 0.5 dex, is consistent 
with the [O/Fe] ratio of about 3 times solar measured in the 
ejecta of the Type II SN1987A (see, e.g., Arnett et al. 1989). 
Thin-disk stars, on the other hand, formed from a gas in which 
a mixture of material from low mass and massive stars was 
present. The steep decline in the thin disk [O/Fe] ratio indicates 
that the chemical enrichment of the ISM was due to both SN la 
and SN II, as proposed originally by Tinsley (1980), and demon- 
strated here, in Fig. [TT] by the satisfactory agreement between 
the data and the Galactic Chemical Evolution (GCE) model by 
Nomoto et al. (2006). 

Inspection of Figs. [8] [9] and QT] reveals that the thick disk 
[O/Fe] ratios decrease from 0.5 to 0.4 with [Fe/H] from [Fe/H] = 
-0.8 to -0.3 and that no knee in the [O/Fe] vs. [Fe/H] trend is 
present at higher metallicities. The small slope can be explained 
by a metallicity dependent yield from SN II because oxygen 
production is favored at low [Fe/H]. Using the nucleosynthe- 
sis yields by Kobayashi et al. (2006, hereafter K06) and a simple 
GCE model (see Appendix|B]), we were able to satisfactorily re- 
produce the slope in the [O/Fe] vs. [Fe/H] relation of the thick 
disk (see Fig.[TTI). Only massive stars were allowed to enrich the 
ISM in our model so, given the good agreement with the data, 
we conclude that the SN II to SN la ratio during the formation 
of the thick disk was large enough to make the contribution of 
SN la to the chemical enrichment of the ISM negligible. This can 
be achieved, for example, by a continuous intense star formation 
during the whole process of thick disk formation. Note that the 
lack of a knee does not exclude SN la events and their contribu- 
tion, albeit small, of metals to the ISM during the formation of 
the thick-disk. 

The [O/Fe] vs. [Fe/H] trend in our simple GCE model is al- 
most insensitive to the choice of star formation rate (e.g., con- 
tinuous, truncated, or exponentially decaying) J3 The metallicity 
dependence of the oxygen yields by K06 lead to a good repre- 
sentation of our observed data. Note, however, that the predicted 
metallicity dependence of the yields for other elements might be 
incompatible with the observational data by, e.g., R06. In par- 
ticular, note the case of Na and Al, for which the yields are 
predicted to be very sensitive to the metallicity. The IMF (ini- 
tial mass function) weighted yields, in the [X/Fe] form, increase 
by about 0.5 and 0.4 dex, respectively, from [Fe/H] ~ -1 to 
[Fe/H] ~ (Fig. 5 in K06), yet the data shows no such depen- 
dence of [Na/Fe] or [Al/Fe] with [Fe/H] in the thick disk (see 
Figs. 15 and 16 in R06). Also, in the case of Mn, the observa- 
tions (Fig. 16 in R06) suggest a strong dependence with metal- 
licity such that the [Mn/Fe] ratio increases by about 0.4 dex from 
[Fe/H] = -1 to solar. Our GCE model for Mn predicted an in- 
crease of only 0.15 dex in this range (see also Fig. 5 in K06). 
Oxygen is probably an exceptionally well modeled element be- 
cause, as stated by K06, it is one of the best described elements 
in nucleosynthesis. 

From Fig.|£l which separates stars with high kinematic prob- 
ability of being thin- or thick-disk stars from those with interme- 
diate kinematics, we conclude that the thin/thick disk separation 
is real, given that including stars with intermediate probabilities 



Partially because we did not include SN la. 
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does not destroy the bimodal distribution. In other words, stars 
with kinematics intermediate between that of the thin and thick 
disks do not belong to another population but are only part of 
the tails of the kinematic distributions, which clearly overlap. 
The discussion in Sect. 16.21 and Fig. [10] allows us to reach the 
same conclusion with more confidence. 

Our results are in qualitative agreement with those by B04, 
who used the forbidden oxygen lines at 630 nm and 636 nm 
in a selected sample of thin- and thick-disk stars. B04, how- 
ever, claim to have found thick disk stars in the range -0.4 < 
[Fe/H] < 0.2 whose [O/Fe] decrease towards the solar value. 
This knee would be the signature of SN la in the abundance pat- 
tern of thick disk stars (see also Feltzing et al. 2003) which we, 
however, cannot confirm. The fact that we do not see the knee is 
not due to the characteristics of our sample (see Fig.[T2land the 
discussion in Sect. [6J. This discrepancy needs to be investigated 
in more detail. 

In view of the age evidence that suggests that thick-disk stars 
formed before thin-disk stars, it is important to note that if a hia- 
tus of star formation occurred between the formation of the two 
disks, it would lead to clearly distinct abundance patterns for 
thin- and thick-disk stars. This, and some of the other arguments 
presented in this discussion, were already proposed in the early 
papers by Gratton et al. (1996) and Fuhrmann (1998). During the 
period of low or no star formation, no more contribution of SN II 
to the ISM is present, leaving AGB stars and possibly SN la (if 
the duration of the hiatus is long enough) as main contributors 
of metals to the ISM. This may explain the sudden decrease in 
[O/Fe] from the thick to the thin disk phase. It does not, how- 
ever, explain the overlap in the [Fe/H] distributions. Additional 
information about the formation history of the Galactic disk is 
thus required. 

Thick disk formation scenarios can be divided into top-down 
and bottom-up models (e.g., Majewski 1993). In the former type 
of models, the thick disk is an evolutionary stage of the present 
thin disk, i.e. an intermediate phase between halo and thin disk 
formation, produced essentially by monolithic collapse, rotation 
and conservation of the angular momentum. Thick-disk stars are 
expected to have an age distribution that connects smoothly with 
that of thin disk stars without an age gap and an abundance gra- 
dient with distance from the galactic plane is predicted. Also, 
in most top-down models, the thin and thick disk abundance 
patterns are expected to be continuous. This is not consistent 
with the observational results, including ours. We cannot probe 
into the existence of a metallicity gradient with distance from 
the Galactic plane, but that has been rejected by Gilmore et al. 
(1995) and Allende Prieto et al. (2006). 

Some of the bottom-up models, on the other hand, start with 
an original thin disk that is later heated by different mechanisms. 
The most accepted scenario in this type of models is heating by 
satellite accretion. Quinn et al. (1993) have shown that the accre- 
tion of a satellite galaxy into a disk galaxy results in a heated disk 
which contains stars from both the original thin disk and the ac- 
creted satellite. They have also shown that if more than one satel- 
lite is accreted, there will be additional heating but only the first 
accretion event is necessary for thick disk formation. Their sim- 
ulation does not include hydrodynamical effects so it is not clear 
what happens to the gas during and after the merger. However, in 
order to explain the observations, it has been suggested that after 
the merger, star formation stops and the gas settles down into a 
new thin disk which is now mixed with the gas from the satellite. 
If the latter is more metal-poor, we expect the new thin disk to 
start with an overall metallicity lower than the maximum metal- 
licity reached by the thick-disk stars, thus explaining the overlap 



in the [Fe/H] distributions of thin- and thick-disk stars. In this 
scenario, due to the hiatus of star formation between thick and 
thin disk formation, the abundance patterns for thin- and thick- 
disk stars are expected to be discrete, as explained above. An 
alternative bottom-up scenario is one in which thick-disk stars 
are formed in pregalactic fragments before or while they merge 
and deposit gas and stars into a disk at large scale heights (e.g., 
Abadi et al. 2003, Brook et al. 2005). After the merger(s), the 
gas settles and cools down, eventually leading to the formation 
of thin-disk stars. 

The details of the different merger scenarios for the for- 
mation of the thick-disk have been summarized by R06 (their 
Sect. 7). In particular, here we emphasize the importance of the 
recent studies by Brook et al. (2004, 2005), who use jV-body 
simulations that take into account, among other things, gas dy- 
namics, star formation, and metal enrichment due to supernovae 
and intermediate mass stars, and reasonably reproduce the obser- 
vational results of R06 as well as ours. Brook et al.'s models pre- 
dict the thin/thick disk separation in the [O/Fe] vs. [Fe/H] plane 
(in fact in the more generic [c/Fe] vs. [Fe/H] plane), extending 
the thick disk [Fe/H] distribution up to values between-0.1 and 
-0.5 depending on the choice of model. In any case, the thick 
disk abundance trend does not connect with that of the thin disk 
at high metallicity, i.e., these models do not predict a knee, in 
good agreement with our results and those by R06. 

Interestingly, additional Galactic and extragalactic evi- 
dence provides even further support to the merger(s) scenario 
(see Freeman & Bland-Hawthorn 2002 for a recent review). 
Abundance patterns for other elements also show, with some 
exceptions, discrete trends for thin- and thick-disk stars (e.g., 
Bensby et al. 2005, R06). Thick disk formation cannot be a nat- 
ural evolutionary phase of non-interacting galaxies, given that 
some disk galaxies seem not to have thick disks (Fry et al. 1999) 
and thick disks have been shown to be more common in interact- 
ing/merging galaxies (Schwarzkopf & Dettmar 2000). Yoachim 
& Dalcanton (2005) have found a galaxy in which the thick disk 
is rotating in the opposite direction with respect to the thin disk, 
a situation that cannot be explained by monolithic collapse. In 
a recent study of thick disks in edge-on galaxies, Yoachim & 
Dalcanton (2006) conclude that their observations can be ex- 
plained by models in which stars form in galactic subunits before 
merging to form a thick disk. Finally, Wyse et al. (2006) present 
an interesting case for streaming stars in the direction of dwarf 
spheroidal galaxies that could be remnants of the past mergers 
that led to the formation of the Galactic thick disk. 

In summary, the abundance patterns derived in this and pre- 
viously published works strongly favor a discrete separation be- 
tween thin- and thick-disk stars. We have demonstrated that an 
intermediate population with properties between those of the 
thin and thick disks does not exist. The systematic differences 
found in the oxygen abundance trends are explained in terms of 
ISM enrichment by SN la and SN II, taking into account the im- 
portance of each of them at different epochs. In particular, we 
showed that the oxygen abundances of thick-disk stars are ex- 
plained by the metallicity dependence of the yields due to mas- 
sive stars only. Stellar dating indicates that the thick disk was 
formed prior to the thin disk, while the abundances suggest a 
hiatus of star formation between thick- and thin-disk formation. 
With the abundance and age evidence, the observed kinematics 
of thin/thick disk stars, and observations of thick disks in other 
galaxies, the merger(s) scenario currently seems to be the most 
reliable model for the formation of the thick disk. It accounts 
for a hiatus of star formation between the formation of the two 
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disks and the abundance and kinematic patterns. Finally, the pre- 
viously proposed existence of a knee in the oxygen abundance 
pattern of the thick disk has been challenged by our results. We 
provide some clues to the reason for the discrepancy but a fuller 
exploration is required, in particular a homogeneous determina- 
tion of oxygen abundances using different indicators and taking 
into account all possible sources of systematic error. 
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Appendix A: Systematic errors in our stellar age 
determination 

In R03, AP04, and this work, the stellar ages were determined 
from the observed stellar parameters and the corresponding val- 
ues predicted by the theoretical isochrones of Bertelli et al. 
(1994). The isochrone points contain the additional information 
required, namely, stellar ages, and also other important quanti- 
ties such as the radius and mass. Thus, a comparison of theoreti- 
cal and observed data provides an estimative of these extra quan- 
tities. The actual procedure involves computing probabilities of 
a given star to have certain ages within the ranges allowed by the 
observed parameters. The probability distribution is then used to 
estimate the age and its error (see Appendix A in AP04a). Here 
we comment on possible systematic errors associated with this 
procedure of age determination. 

Given the impact of a-element enhancement on stellar evo- 
lution calculations (e.g., Kim et al. 2002), the use of models 
with solar-scaled chemical compositions (like those by Bertelli 
et al. 1994) introduces systematic errors in the age determina- 
tion from isochrones. We compared the ages determined with 
and without this effect. In addition, since the age determination 
from isochrones is subject to statistical biases (e.g., Pont & Eyer 



2004, J0rgensen & Lindegren 2005), we compared our ages with 
those derived using the Pont & Eyer (2004) method, which uses 
Bayesian statistics and a procedure that attempts to eliminate 
biases in standard isochrone age determinations. Note also that 
the Pont & Eyer method uses the synthetic color-magnitude di- 
agrams computed by Aparicio & Gallart (2004), which, in addi- 
tion to the Bertelli et al. results, is also based on isochrones that 
consider a-element enhancement. 

Using the empirical prescription by DeglTnnocenti et al. 
(2005), in which the a-element enhancement effects on the 
isochrones are mimicked by increasing the [Fe/H] values, as 
originally proposed by Salaris et al. (1993)0 we find that the 
ages obtained without a-element enhancement are underesti- 
mated with respect to those obtained with it. The differences be- 
tween ages computed with and without the a-element enhance- 
ment effects are more important for the oldest stars. The differ- 
ences are as large as 4 Gyr for a 13 Gyr star and 2 Gyr for an 
8 Gyr star. 

On the other hand, our isochrone ages are in satisfactory 
agreement with those obtained from the Pont & Eyer (2004) 
method, considering the uncertainties of both age estimates. On 
average, below 6 Gyr, a constant difference of about 1 Gyr is 
observed, with our mean ages being smaller. From about 10 
to 15 Gyr, the differences with the Pont & Eyer ages increase 
with age, but this time in the opposite direction (our adopted 
ages are larger). The largest differences are for the oldest stars, 
whose ages may be overestimated by about 3 Gyr when using 
our method. Even though this may be related to the 'terminal 
age bias' discussed by Pont & Eyer (2004), it is more likely to 
be explained by systematic differences in the isochrones (note, 
again, that the input data used by the Pont & Eyer method in- 
cludes isochrones computed with a-element enhancement). 

Tentatively, we conclude that using the Pont & Eyer (2004) 
method and/or considering the a-element enhancement effect 
would result in ages substantially lower than those obtained from 
our method for the oldest stars. Although these studies are aim- 
ing at a substantial improvement to the methods of stellar dat- 
ing, they need to be tested in more detail. Note, for example, 
that the a-element enhancement leads to a thick disk younger 
than 10 Gyr. Although this is in disagreement with other re- 
cent studies of the thick disk (but note that Allende Prieto et al. 
2006 conclude that thick-disk stars as young as 8 Gyr exist), it 
solves the problem that, according to our calculations, a few (4) 
of our stars (and a similar number in the Reddy et al. studies) 
result with lower limits in their ages that are greater than the 
upper limit on stellar ages (the age of the universe) imposed by 
the Wilkinson Microwave Anisotropy Probe (WMAP) measure- 
ments (13.7 ± 0.2 Gyr, Krauss 2003). 

It should be noted that these effects have a smaller impact 
on our log g determination. For example, if we adopt a-element 
enhanced isochrones our log g values for dwarf stars, on average, 
increase only by about 0.05 dex, thus introducing a very small 
systematic error in the abundances. 

Appendix B: A simple GCE model for the thick disk 

Simple homogeneous Galactic Chemical Evolution (GCE) mod- 
els have been described in detail by, e.g., Tinsley (1980), Clayton 



17 Here the original (percentile) metallicity Z is increased to Z = 
Zo(af„ + b), where /„ = 10 [ff/Fel and a, b are coefficients that depend on 
the adopted solar abundances. The a, b values we used, which are given 
by DeglTnnocenti et al. (2005), correspond to the solar abundances by 
Asplund et al. (2005). 
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t (Gyr) 


[Fe/H] 


SFR 
[V/te\ 


= const. 
|Mn/reJ 


oc e 


-1/lGyr 
|Mn/reJ 


oc e 
[0/te\ 


-f/4Liyr 

[Mn/reJ 


Truncated at 1 Gyr 
[O/Fe] [Mn/Fe] 


0.05 


-1.66 


0.601 


-0.536 


0.601 


-0.536 


0.601 


-0.536 


0.601 


-0.536 


0.08 


-1.52 


0.582 


-0.539 


0.582 


-0.539 


0.582 


-0.539 


0.582 


-0.539 


0.13 


-1.38 


0.564 


-0.542 


0.564 


-0.542 


0.564 


-0.542 


0.564 


-0.542 


0.20 


-1.24 


0.544 


-0.545 


0.545 


-0.544 


0.545 


-0.545 


0.544 


-0.545 


0.32 


-1.10 


0.523 


-0.542 


0.525 


-0.542 


0.523 


-0.542 


0.523 


-0.542 


U.JU 




U.470 






U.JJU 




U.J J J 


0.498 


-0.535 


0.79 


-0.82 


0.473 


-0.518 


0.482 


-0.523 


0.475 


-0.519 


0.473 


-0.518 


1.26 


-0.68 


0.445 


-0.489 


0.462 


-0.504 


0.450 


-0.492 


0.458 


-0.504 


2.00 


-0.54 


0.427 


-0.455 


0.450 


-0.485 


0.433 


-0.463 


0.458 


-0.505 


3.16 


-0.40 


0.415 


-0.417 


0.442 


-0.470 


0.421 


-0.431 


0.454 


-0.502 


5.01 


-0.26 


0.406 


-0.378 


0.435 


-0.462 


0.413 


-0.402 


0.455 


-0.508 



Table B.l. Results of our thick disk GCE model for different choices of the SFR. 



(1984), Rana (1991), Pagel (1997), and Matteucci (2001). Here 
we present a very simple GCE model in which several assump- 
tions are made to simplify the problem but remains a reason- 
ably good picture to study the general properties of the chemical 
evolution of the thick disk. Readers should be warned that this 
model does not follow the evolution of the gas fraction and does 
not take into account infall (hence it is not suitable to study the 
thin/thick disk transition). Our goal with this calculation is to 
show what the effect of a metallicity dependent yield of SN II 
on the evolution of abundance ratios in an isolated interstellar 
medium (ISM) is. Although it can be extended to include SN la 
yields, we computed the case in which only massive stars con- 
tribute to the enrichment of the ISM. Therefore, this model can- 
not be used to study the G-dwarf problem either. More elabo- 
rated models for the thin/thick disk chemical evolution can be 
found elsewhere (e.g., Chiappini et al. 1997, Brook et al. 2005). 

We computed the number density of atoms of an element X 
in the ISM at a time t from the following formula: 



N x (t) 



Jo Jt' 



^(f'OOtMtf', t")]Y x [M(t', t"),Z(t")]dt"dt' 



where *P is the star formation rate (SFR), <t> the initial mass func- 
tion (IMF), Yx the yield of X due to supernovae, M the stellar 
mass and Z the metallicity of the ISM. The limits of integration 
and f max are related to the minimum (M mm ) and maximum 
(M max ) masses of the stars that enrich the ISM in the model: 



Cn = max(f'-MCf,0), and 



to - 1,2,4 Gyr, and a SFR truncated at 1 Gyr but the final results 
were very similar when compared to the observed data. Since in 
this model we have no control on the total amount of gas we 
had to assume an age-metallicity relation in order to compute 
abundance ratios. The following age-metallicity relation: 



Z(f) 



3.56 x lO^f - 7 



which is a reasonably good approximation to the theoretical 
results by Chiappini et al. (1997), and also to our observed 
data (see Fig. [T3l ). was adopted. We then assumed [Fe/H] = 
log(Z/Z Q ). Finally, the abundance ratios were computed as fol- 
lows: 



[X/Y] 



& \N Y (t)J 



(A x - A F ) . 



where the difference (Ax - Ay) Q was obtained from the so- 
lar abundances. Note that by taking the abundance ratio we 
eliminate the proportionality constants in the IMF and SFR, 
which further simplifies the problem. Note also that if the age- 
metallicity relation is inaccurate, only the time dependence of 
the ratios would be affected, but not the abundance trends, e.g., 
the [O/Fe] vs. [Fe/H] relations, as long as the SFR is a slowly 
varying function of time (in particular for the case of a constant 
SFR). 

The results and implications of this model for the thick disk 
are discussed in Sect. [8] and are shown in Fig. QT|for the case 
of constant SFR. The [O/Fe] and [Mn/Fe] results are given in 
Table iBll for different choices of the SFR. 



4ax = f ~ IqM, 



-2.5 
max 



where Iq comes from the relation adopted for the lifetime r of a 
star, t = l()M~ 2 5 , that results from assuming the luminosity of 
massive stars to be roughly proportional to M 3 5 (e.g., Demircan 
& Kahraman 1991). Clearly, lo is the lifetime of the Sun (ap- 
proximately 9 Gyr) and M is given in solar masses. The mass of 
a star born at t" and exploding at t' is then 



M(t' \t") = Z°V 



t")- QA 



We used the mass- and metallicity-dependent yields (Yx) by 
K06 (see also Nomoto et al. 2006), who provide yields for nor- 
mal SN II as well as hypernovae. We assumed the ratio of hy- 
pernovae to normal SN II to be 0.5, a value that reproduces very 
well abundance trends in more complex GCE models (K06). 

The Salpeter (1955) IMF, <D(M) oc AT 2 35 , was adopted. 
Then, 0[M(f', f")] oc (f - t") 0M , and the integration is straight- 
forward. For the SFR we used ^(t) = const., ^(f) oc e~' lla , where 
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Table 4. Basic and kinematic data. The full version of this table, which will be available online at the CDS, can be obtained from the authors upon 
request. 



HIP 


HD 


V 


K 


v r 


U 


V 


W 


Pi 


Pi 


f 3 






mag 


mas 


km s -1 


km s _1 


km s _1 


km s~ ] 








171 


224930 


5.75 


80.6 ±3.0 


-36.9 ± 0.4 


-8.0 ± 0.7 


-73.1 ± 1.7 


-31.1 ± 2.0 


0.23 


0.76 


0.01 


394 


225239 


6.11 


27.2 ±4.3 


4.6 ± 0.3 


-124.8 ± 19.4 


-51.3 ±8.7 


-10.9 ± 1.5 


0.27 


0.71 


0.02 


475 


70 


8.22 


21.0± 1.1 


-28.4 ± 0.2 


6.8 ±0.4 


-31.4 ± 0.4 


-36.7 ± 2.0 


0.93 


0.07 


0.00 


493 


101 


7.46 


26.2 ± 0.8 


-45.6 ±0.1 


45.8 ± 1.1 


-34.2 ±0.1 


17.0 ±0.5 


0.95 


0.05 


0.00 


522 


142 


5.71 


39.0 ± 0.6 


5.3 ±0.3 


-57.4 ± 1.0 


-37.5 ± 0.6 


-15.3 ±0.3 


0.93 


0.07 


0.00 


530 


153 


8.36 


8.1 ±0.9 


-31.8 ±0.3 


-28.2 ± 4.5 


-48.5 ± 2.4 


2.8 ± 0.9 


0.92 


0.08 


0.00 
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Table 5. Line data. The Van der Waals broadening cross section (<r) for an atom-perturber relative velocity vo = 10 6 cm s -1 and velocity parameter 
a are given. The a values are given in atomic units (1 a.u. = 2.80 x 1(T 17 cm 2 ) The damping constants j(, (FWHM) per H-atom at T = 10 4 K 
and C 6 , which are used as input data by some spectrum synthesis codes, can be computed from the values given in this table using the following 
relations (see, e.g., Gray 1992, Anstee & O'Mara 1995): 

t = 2(if 2 r(^)v.( t )- = i7v3/ 5 cf , 

where v = (SkT/nj/) [l2 , ji 1 = (m x l + m H [ ), m x is the mass of an atom of the element X, k the Boltzmann constant, and T the gamma function. For 
the lines marked with a t the damping constants correspond to the modified Unsold approximation (C 6 multiplied by 2). The last column gives the 
solar equivalent widths. 



Wavelength 


Species 


EP 


logg/ 


cr 


a 


EW 


A 




eV 




a.u. 




mA 


4620.52 


Fen 


2.83 


-3.21 


185 


0.306 


55.6 


4629.34 


Fe 11 


2.81 


-2.28 


178 


0.257 


98.5 


4630.12 


Fei 


2.28 


-2.52 


416 


0.254 


73.6 


4635.85 


Fei 


2.85 


-2.34 


277 


0.252 


53.8 


4683.56 


Fei 


2.83 


-2.41 


218 


0.263 


56.2 


4690.14 


Fei 


3.69 


-1.61 


792 


0.283 


54.3 


4745.80 


Fei 


3.65 


-1.27 


621 


0.400 t 


78.0 


4749.95 


Fei 


4.56 


-1.24 


634 


0.275 


35.3 


4779.44 


Fei 


3.42 


-2.16 


317 


0.261 


40.6 


4788.76 


Fei 


3.24 


-1.73 


238 


0.249 


66.1 


4799.41 


Fei 


3.64 


-2.13 


353 


0.241 


33.9 


4808.15 


Fei 


3.25 


-2.69 


297 


0.274 


27.8 


4961.91 


Fei 


3.63 


-2.19 


375 


0.256 


25.9 


4962.57 


Fei 


4.18 


-1.18 


804 


0.400 t 


55.3 


5044.21 


Fei 


2.85 


-2.03 


713 


0.238 


73.5 


5054.64 


Fei 


3.64 


-1.98 


353 


0.313 


38.3 


5145.09 


Fei 


2.20 


-3.08 


351 


0.272 


52.6 


5187.91 


Fei 


4.14 


-1.26 


724 


0.400 t 


63.4 


5197.58 


Fen 


3.23 


-2.22 


180 


0.247 


81.4 


5217.92 


Fei 


3.64 


-1.72 


542 


0.400 t 


50.7 


5228.38 


Fei 


4.22 


-1.19 


809 


0.278 


61.2 


5234.63 


Fe 11 


3.22 


-2.18 


180 


0.249 


82.3 


5253.46 


Fei 


3.28 


-1.57 


849 


0.229 


81.9 


5264.81 


Fe 11 


3.23 


-3.13 


186 


0.300 


40.5 


5285.13 


Fei 


4.43 


-1.54 


1046 


0.282 


30.6 


5288.52 


Fei 


3.69 


-1.51 


353 


0.297 


60.0 


5293.96 


Fei 


4.14 


-1.77 


701 


0.400 t 


31.5 


5295.31 


Fei 


4.42 


-1.59 


1014 


0.281 


28.1 


5322.04 


Fei 


2.28 


-2.89 


341 


0.236 


60.8 


5373.71 


Fei 


4.47 


-0.74 


1044 


0.282 


63.8 


5379.57 


Fei 


3.69 


-1.51 


363 


0.249 


61.8 


5386.33 


Fei 


4.15 


-1.67 


930 


0.278 


32.7 


5414.07 


Fe 11 


3.22 


-3.58 




0.303 


28.8 


5432.95 


Fei 


4.45 


-0.94 


974 


0.280 


70.9 


5441.34 


Fei 


4.31 


-1.63 


807 


0.278 


36.8 


5464.28 


Fei 


4.14 


-1.58 


668 


0.400 t 


38.1 


5472.71 


Fei 


4.21 


-1.52 


752 


0.215 


45.8 


5473.90 


Fei 


4.15 


-0.72 


738 


0.241 


84.5 


5483.10 


Fei 


4.15 


-1.45 


737 


0.241 


47.9 


5487.15 


Fei 


4.42 


-1.43 


908 


0.279 


41.0 


5522.45 


Fei 


4.21 


-1.45 


744 


0.215 


43.4 


5525.12 


Fe 11 


3.27 


-3.94 


184 


0.280 


15.1 


5525.54 


Fei 


4.23 


-1.12 


749 


0.238 


55.8 


5543.94 


Fei 


4.22 


-1.04 


742 


0.238 


63.5 


5546.51 


Fei 


4.37 


-1.21 


825 


0.278 


52.9 


5560.21 


Fei 


4.43 


-1.09 


895 


0.278 


51.9 


5584.76 


Fei 


3.57 


-2.22 


296 


0.271 


37.9 


5618.63 


Fei 


4.21 


-1.27 


732 


0.214 


50.5 


5619.60 


Fei 


4.39 


-1.60 


808 


0.277 


35.7 


5633.95 


Fei 


4.99 


-0.23 


635 


0.270 


69.9 


5635.82 


Fei 


4.26 


-1.79 


928 


0.279 


34.0 


5638.26 


Fei 


4.22 


-0.77 


730 


0.235 


78.2 


5650.71 


Fei 


5.09 


-0.86 


1457 


0.400 t 


42.0 


5653.87 


Fei 


4.39 


-1.54 


792 


0.277 


38.9 


5679.02 


Fei 


4.65 


-0.82 


1106 


0.291 


60.0 


5701.54 


Fei 


2.56 


-2.22 


361 


0.237 


84.8 


5705.46 


Fei 


4.30 


-1.40 


744 


0.231 


38.8 


5717.83 


Fei 


4.28 


-1.03 


758 


0.209 


66.6 


5731.76 


Fei 


4.26 


-1.20 


727 


0.232 


59.4 


5775.08 


Fei 


4.22 


-1.30 


648 


0.400 t 


59.7 


5793.92 


Fei 


4.22 


-1.60 


714 


0.231 


35.0 


5806.73 


Fei 


4.61 


-0.95 


985 


0.281 


55.0 


5809.22 


Fei 


3.88 


-1.74 


956 


0.244 


52.0 


5852.22 


Fei 


4.55 


-1.23 


896 


0.276 


41.6 


5856.09 


Fei 


4.29 


-1.46 


404 


0.264 


33.8 


5905.67 


Fei 


4.65 


-0.69 


994 


0.282 


60.0 


5909.97 


Fei 


3.21 


-2.60 


716 


0.243 


40.7 


5916.25 


Fei 


2.45 


-2.99 


341 


0.238 


57.0 


5927.79 


Fei 


4.65 


-0.99 


984 


0.281 


42.7 


5929.68 


Fei 


4.55 


-1.31 


864 


0.275 


40.6 


5930.18 


Fei 


4.65 


-0.17 


983 


0.281 


90.1 


5934.65 


Fei 


3.93 


-1.07 


959 


0.247 


79.7 


5952.72 


Fei 


3.98 


-1.34 


999 


0.252 


67.6 
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Table 5. continued. 



Wavelength 


Species 


EP 


logs/ 


cr 


a 


EW 


A 




eV 




a.u. 




mA 


6003.01 


Fei 


3.88 


-1.06 


898 


0.241 


87.9 


6027.05 


Fei 


4.08 


-1.09 


567 


0.400 t 


65.6 


6056.00 


Fei 


4.73 


-0.40 


1029 


0.286 


74.8 


6079.01 


Fei 


4.65 


-1.02 


920 


0.276 


47.9 


6082.71 


Fei 


2.22 


-3.57 


306 


0.271 


36.2 


6085.26 


Fei 


2.76 


-3.05 


344 


0.260 


44.0 


6093.64 


Fei 


4.61 


-1.40 


866 


0.274 


31.7 


6096.67 


Fei 


3.98 


-1.83 


963 


0.250 


39.0 


6127.91 


Fei 


4.14 


-1.40 


574 


0.400 t 


51.0 


6151.62 


Fei 


2.18 


-3.30 


277 


0.263 


50.9 


6165.36 


Fei 


4.14 


-1.46 


570 


0.400 f 


46.4 


6170.51 


Fei 


4.80 


-0.38 


1057 


0.289 


86.1 


6173.34 


Fei 


2.22 


-2.88 


281 


0.266 


69.7 


6180.20 


Fei 


2.73 


-2.63 


335 


0.259 


58.4 


6187.99 


Fei 


3.94 


-1.62 


903 


0.244 


48.9 


6200.31 


Fei 


2.61 


-2.44 


350 


0.235 


75.7 


6213.43 


Fei 


2.22 


-2.48 


280 


0.265 


84.7 


6229.23 


Fei 


2.85 


-2.83 


350 


0.250 


39.4 


6232.64 


Fei 


3.65 


-1.22 


446 


0.400 t 


87.1 


6265.13 


Fei 


2.18 


-2.55 


274 


0.261 


88.6 


6270.23 


Fei 


2.86 


-2.54 


350 


0.249 


54.5 


6297.79 


Fei 


2.22 


-2.74 


278 


0.264 


76.5 


6322.69 


Fei 


2.59 


-2.43 


345 


0.238 


78.1 


6344.15 


Fei 


2.43 


-2.92 


327 


0.246 


66.6 


6380.74 


Fei 


4.19 


-1.32 


562 


0.400 t 


54.3 


6419.95 


Fei 


4.73 


-0.24 


885 


0.274 


87.6 


6432.68 


Fen 


2.89 


-3.57 


169 


0.204 


43.9 


6481.87 


Fei 


2.28 


-2.98 


308 


0.243 


66.6 


6518.37 


Fei 


2.83 


-2.52 


336 


0.251 


59.9 


6593.87 


Fei 


2.43 


-2.42 


321 


0.247 


88.0 


6597.56 


Fei 


4.80 


-0.97 


893 


0.276 


44.7 


6609.11 


Fei 


2.56 


-2.69 


335 


0.245 


68.6 


6703.57 


Fei 


2.76 


-3.06 


320 


0.264 


38.5 


6750.15 


Fei 


2.42 


-2.62 


335 


0.241 


75.8 


6752.71 


Fei 


4.64 


-1.22 


778 


0.274 


37.6 


6810.26 


Fei 


4.61 


-1.00 


873 


0.275 


50.4 


6828.59 


Fei 


4.64 


-0.82 


768 


0.277 


59.1 


6841.34 


Fei 


4.61 


-0.71 


759 


0.267 


71.9 


6842.69 


Fei 


4.64 


-1.22 


896 


0.279 


40.7 


6843.66 


Fei 


4.55 


-0.83 


736 


0.216 


61.7 


6855.16 


Fei 


4.56 


-0.74 


637 


0.400 t 


76.7 


6858.15 


Fei 


4.61 


-0.94 


765 


0.211 


52.8 


6945.21 


Fei 


2.42 


-2.48 


331 


0.243 


84.3 


6978.85 


Fei 


2.48 


-2.50 


337 


0.241 


81.6 


6999.88 


Fei 


4.10 


-1.46 


845 


0.244 


58.2 


7022.95 


Fei 


4.19 


-1.15 


912 


0.245 


68.4 


7038.22 


Fei 


4.22 


-1.20 


933 


0.247 


66.0 


7090.38 


Fei 


4.23 


-1.11 


934 


0.248 


70.8 


7132.99 


Fei 


4.08 


-1.65 


473 


0.400 t 


44.5 


7222.39 


Fe ii 


3.89 


-3.26 


187 


0.301 


16.3 


7224.49 


Fe ii 


3.89 


-3.20 


187 


0.301 


22.7 


7401.69 


Fei 


4.19 


-1.55 


481 


0.400 t 


41.6 


7418.67 


Fei 


4.14 


-1.38 


469 


0.400 t 


48.8 


7449.33 


Fen 


3.89 


-3.27 


186 


0.301 


23.5 


7515.83 


Fe ii 


3.90 


-3.39 


186 


0.301 


13.6 


7583.79 


Fei 


3.02 


-1.88 


367 


0.238 


85.7 


7710.36 


Fei 


4.22 


-1.11 


470 


0.400 t 


68.6 


7711.72 


Fen 


3.90 


-2.50 


185 


0.300 


45.3 


7723.21 


Fei 


2.28 


-3.62 


304 


0.242 


44.8 


7771.94 


Oi 


9.15 


0.37 


453 


0.234 


70.6 


7774.17 


Oi 


9.15 


0.22 


453 


0.234 


60.6 


7775.39 


Oi 


9.15 


0.00 


453 


0.234 


48.8 
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Table 6. Fundamental parameters, abundances, and ages. The full version of this table, which will be available online at the CDS, can be obtained 
from the authors upon request. 



HIP 


HD 


Teg 


logg 


V, 


[Fe/H] F( . , 


[Fe/H] Fc n 


[Fe/H] 


[0/H] LTE 


[0/H] NLXE 


log Age 






K 


cm s~ 2 


km s _1 


dex 


dex 


dex 


dex 


dex 


yr 


171 


224930 


5405 


4.60 


0.86 


-0.86 ± 0.06 


-0.76 ± 0.06 




-0.10 ±0.01 


-0.16 ±0.03 




394 


225239 


5537 


3.74 


1.16 


-0.56 ± 0.04 


-0.52 ± 0.06 


-0.51 ±0.04 


-0.05 ± 0.03 


-0.23 ± 0.01 




475 


70 


5660 


4.47 


0.84 


-0.43 ± 0.04 


-0.33 ± 0.04 


-0.35 ± 0.03 


-0.02 ± 0.01 


-0.08 ± 0.02 




493 


101 


5858 


4.45 


1.01 


-0.28 ± 0.04 


-0.21 ±0.03 


-0.21 ± 0.03 


-0.04 ± 0.08 


-0.11 ±0.05 




522 


142 


6177 


4.25 


1.72 


-0.05 ± 0.08 


-0.06 ± 0.09 


-0.02 ± 0.06 


0.27 ± 0.06 


0.12 ±0.04 





530 153 5779 3.80 1.35 -0.11 ±0.06 -0.09 ± 0.04 -0.08 ± 0.04 0.34 ± 0.04 0.14 ±0.02 9.550 



